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Background: Type la supernovae contribute significantly to the nucleosynthesis of many Fe-group and 
intermediate-mass elements. However, the robustness of nucleosynthesis obtained via models of this class of 
explosions has not been studied in depth until now. 

Purpose: We explore the sensitivity of the nucleosynthesis resulting from thermonuclear explosions of massive 
white dwarfs with respect to uncertainties in nuclear reaction rates. We lay particular emphasis on indentifying 
the individual reactions rates that most strongly affect the isotopic products of these supernovae. 

Method: We have adopted a standard one-dimensional delayed detonation model of the explosion of a 
Chandrasekhar-mass white dwarf, and have post-processed the thermodynamic trajectories of every mass-shell 
with a nucleosynthetic code in order to obtain the chemical composition of the ejected matter. We have con- 
sidered increases (decreases) by a factor of ten on the rates of 1196 nuclear reactions (simultaneously with their 
inverse reactions) repeating the nucleosynthesis calculations after modification of each reaction rate pair. We 
have computed as well hydrodynamic models for different rates of the fusion reactions of 12 C and of le O. From 
the calculations we have selected the reactions that have the largest impact on the supernova yields, and we 
have computed again the nucleosynthesis using two or three alternative prescriptions for their rates, taken from 
the JINA REACLIB database. For the three reactions with the largest sensitivity we have analyzed as well the 
temperature ranges where a modifications of their rates has the strongest effect on nucleosynthesis. 

Results: The nucleosynthesis resulting from the Type la supernova models is quite robust with respect to vari- 
ations of nuclear reaction rates, with the exception of the reaction of fusion of two 12 C nuclei. The energy of the 
explosion changes by less than ~ 4% when the rates of the reactions 12 C + 12 C or 16 + ie O are multiplied by 
a factor of xlO or xO.l. The changes in the nucleosynthesis due to the modification of the rates of these fusion 
reactions are as well quite modest, for instance no species with a mass fraction larger than 0.02 experiences a 
variation of its yield larger than a factor of two. We provide the sensitivity of the yields of the most abundant 
species with respect to the rates of the most intense reactions with protons, neutrons, and alphas. In general, 
the yields of Fe-group nuclei are more robust than the yields of intermediate-mass elements. Among the species 
with yields larger than 10 -8 Mq, 35 S has the largest sensitivity to the nuclear reaction rates. It is remark- 
able that the reactions involving elements with Z > 22 have a tiny influence on the supernova nucleosynthesis. 
Among the charged particle reactions, the most influential on supernova nucleosynthesis are 30 Si + pfi 31 P + 7, 
20 Ne + a «=t 24 Mg + 7, and 24 Mg + a <=t 27 A1 + p. The temperatures at which a modification of their rate has a 
larger impact are in the range 2 < T < 4 GK. 

Conclusions: The explosion model (i.e., the assumed conditions and propagation of the flame) chiefly determines 
the element production of Type la supernovae, and derived quantities like their luminosity, while the nuclear 
reaction rates used in the simulations have a small influence on the kinetic energy and final chemical composition 
of the ejecta. Our results show that the uncertainty in individual thermonuclear reaction rates cannot account 
for discrepancies of a factor of two between isotopic ratios in Type la supernovae and those in the solar system, 
especially within the Fe-group 
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I. INTRODUCTION 

Thanks to their high luminosity, Type la supernovae 
(SNIa) are used routinely as standard candles to mea- 
sure cosmological distances. They are instrumental to 
our current understanding of the Universe, providing ev- 
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idence for its accelerated expansion PHf|. Type la su- 
pernovae play also an important role in the chemical 
evolution of galaxies, being responsible for most of the 
Fe-group elements and smaller amounts of Silicon, Sul- 
fur, Argon, and Calcium (see e.g. @-IU). The elemen- 
tal composition is evident in optical and infrared spectra 
recorded from days to months after the explosion (see 
e.g. (ToL fTl|) and in X-ray spectra of their remnants visi- 
ble for hundreds of years (for a review see [l2j ). Finally, 
SNIa are one of the key targets for 7-ray astronomy, as a 
source of a variety of radioactive isotopes (see e.g. |15|). 
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This is the first paper of a series in which we will study 
the sensitivity of the nucleosynthesis produced in SNIa 
with respect to uncertainties in nuclear data. In this 
paper, we study the sensitivity to variations in rates of 
thermonuclear reactions (fusion reactions, radiative cap- 
tures, and transfer reactions). In forthcoming publica- 
tions, we will study the sensitivity to uncertainties in 
nuclear masses and in weak interaction rates. Studies of 
the effect of nuclear data uncertainties in different astro- 
physical scenarios have been published from time to time 
during the past few decades, e.g. [T3 - l2]| to cite only a 
few, although none of them has dealt with SNIa. These 
works followed different methodologies to test the impact 
of nuclear reaction rates. For instance, ref. [l9| varied the 
rate of individual nuclear reactions relevant for 44 Ti nu- 
cleosynthesis in order to determine which reactions were 
a prime target for the experimental measurement of their 
cross sections. On the other hand ref. [20| designed a 
numerical experiment to measure the uncertainty of the 
nucleosynthesis of nova explosions. To this end, they 
followed a Monte Carlo approach in which they varied 
simultaneously by random factors all the reaction rates 
in their network. The focus of this second approach was 
on the final nova nucleosynthesis rather than in deter- 
mining the individual reactions that are most influential. 
Finally, ref. [2l[ used theoretical nuclear reaction rates 
based on four different nuclear mass models to determine 
their impact on the r-process abundances. In this case, 
the emphasis was on testing different nuclear models. In 
the present work, we wish to determine the individual nu- 
clear reactions most influential on the nucleosynthesis of 
SNIa, hence we will follow the same strategy as ref. . 

At present, the favored model of SNIa is the thermonu- 
clear explosion of a carbon-oxygen white dwarf (WD) 
near the Chandrasekhar mass that accretes matter from a 
companion star in a close binary system [22| ■ Other mod- 
els, such as the sub-Chandrasekhar models or the double 
degenerate scenario, although not completely ruled out, 
either have difficulties in explaining the gross features of 
the spectrum and light curve of normal SNIa, or face 
severe theoretical objections, see e.g. [23T4281 ] . Super- 
Chandrasekhar models have been proposed to explain a 
few overluminous SNIa (29l - [33| but, given the scarcity 
of observations despite of their high intrinsic luminos- 
ity they are thought to represent at most a few percent 
of all SNIa explosions. Moreover, the properties of the 
progenitors of super-Chandrasekhar SNIa and the explo- 
sions themselves are not well understood. Thus, we will 
concentrate our efforts on the study of a reference SNIa 
Chandrasekhar- mass model [13] . 

Even though the hydrostatic evolution of SNIa pro- 
genitors lasts for several Gyrs while the thermonuclear 
explosion lasts for a few seconds at most, the outcome is 
nearly independent of the history of the white dwarf prior 
to its explosive ignition. This fact is commonly denoted 
as 'stellar amnesia' [35| . The only link between the white 
dwarf at ignition time and its previous evolution comes 
through its chemical composition ( 12 C, 16 , 22 Ne, and 



other trace species) and the distribution of hot spots that 
are the seeds of the emerging thermonuclear flame. The 
influence of uncertain reaction rates, specifically that of 
the reaction 12 C(a, 7) 16 0, on the chemical composition 
of massive white dwarfs has been studied by ref. [361 ] , who 
found that the central C/O ratio might vary by a factor 
of ~ 13 at ignition time. On the other hand, the effect of 
different C/O ratios on supernova luminosity and nucle- 
osynthesis was studied in ref. [37[, therefore accounting 
implicitly for a variation on the rate of a capture on 12 C. 
They found that the C/O ratio can have a sizeable im- 
pact on the ejecta composition. On the contrary, ref. [38[ 
reached the opposite conclusion after analyzing the same 
problem with their three-dimensional deflagration mod- 
els of SNIa. 

Prior to the SNIa explosion there is a phase of carbon 
simmering that lasts ~ 1000 yrs and involves tempera- 
tures below 10 9 K. During this phase the neutron excess 
of matter can be raised due to electron captures on 13 N 
and 23 Na. The leading thermonuclear reactions during 
carbon simmering are, aside from 12 C+ 12 C, reactions 
that participate in the transmutation of 12 C into 16 O: 
12 C(p, 7 ) 13 N, 12 C(n, 7 ) 13 C, and 13 C(a,n) 16 0. However, 
the timescale of neutronization is controlled by the 12 C 
fusion reaction and the rate of electron captures [HI, H(| ■ 
Thus, we do not expect that a modification of the rates of 
radiative captures and transfer reactions can affect appre- 
ciably the neutronization of the white dwarf and, hence, 
the final supernova composition. In this work, we will 
consider only modifications of the thermonuclear reac- 
tion rates during the explosive phase of the supernova. 

The temperature range relevant for explosive nucle- 
osynthesis in SNIa is approximately 10 9 K to 10 10 K. 
However, at densities and temperatures in excess of 
~ 10 8 g cm -3 and ~ 5.5 x 10 9 K nuclei attain a nuclear 
statistical equilibrium state (NSE) in which the chemical 
composition, for given temperature, density and electron 
mole fraction, is determined by nuclear bulk properties 
(masses and partition functions), i.e. it does not depend 
on the reaction rates. In these conditions, NSE erases 
any imprint of the previous thermodynamic evolution of 
matter, and reaction rates do not play any role until mat- 
ter leaves NSE (freeze-out process). The minimum tem- 
perature relevant for nucleosynthesis in SNIa depends on 
the type of combustion front. For a detonation, a shock 
heats the fuel to temperatures > 2 x 10 9 K, the pre- 
cise value depending mainly on density, before nuclear 
reactions start modifying the chemical composition. On 
the other hand, the process of combustion within a sub- 
sonic flame presents two different phases. Below a critical 
temperature, T cr j t ~ 2-5 x 10 9 K, the matter tempera- 
ture is set by heat diffusion from the hot ashes, while 
above T cr j t the nuclear energy released by combustion 
dominates over heat diffusion. Thus, we do not expect 
modifying the thermonuclear reaction rates below ~ 1— 
2 x 10 9 K to have an impact on the final chemical com- 
position. 

The plan of the paper is as follows. In the next Section, 
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we detail the methodology used to achieve our goals. We 
describe the post-processing code used to integrate the 
nuclear evolutionary equations, the characteristics of our 
reference SNIa model, the selection of the nuclear reac- 
tions to test for variations in their rates, and the ways in 
which we have modified these rates. In Section IIII| we 
present the results of the sensitivity study with respect 
to the fusion reactions of 12 C, 16 0, and the 3a reaction, 
which are the reactions that rule the initial steps of ther- 
monuclear combustion in SNIa. We test modifications 
of the first two reaction rates for effects on the propaga- 
tion of the flame during a SNIa explosion. In Section ITVl 
we present the results of the sensitivity study with re- 
spect to thermonuclear reaction rates involving protons, 
neutrons, and a particles. We have followed different 
strategies in modifying these reaction rates, using either 
a fixed enhancement factor or a temperature dependent 
one. We have tested as well the use of different prescrip- 
tions for the most influential reaction rates, taken from 
recent literature. For a few reactions we have explored 
the temperature range where a modification of their rates 
have a stronger impact on the supernova yields. Finally, 
in Section [Vj we summarize and give our conclusions. 



II. METHODOLOGY 

A. Integration of the nuclear evolutionary 
equations 

We have computed the chemical composition of a refer- 
ence SNIa model with the nucleosynthetic code CRANK 
(Code for the Resolution of an Adaptive Nuclear net- 
worK). CRANK is a post-processing code that integrates 
the temporal evolution of a nuclear network for given 
thermal and structural (density) time profile, and ini- 
tial composition. We have selected the nuclear reactions 
that contribute most to the synthesis of abundant species. 
Then, we have recomputed the nucleosynthesis modifying 
the rate of each one of the selected reactions. 

The inputs to CRANK are the nuclear data and the 
thermodynamic trajectories, as a function of time, of 
each mass shell of the supernova model. The evolution- 
ary equations for the nuclear composition follow the time 
evolution of the molar fraction, Yi , or abundance of each 
species until the temperature falls below 10 8 K, after 
which time the chemical composition is no longer sub- 
stantially modified. The nuclear network is integrated 
with an implicit, iterative method with adaptive time 
steps. The iterative procedure ends when the molar 
abundances of all species with Yi > 10~ 14 mol g _1 have 
converged to better than a relative variation of 10~ 6 . 

The nuclear species present in the network are dynam- 
ically determined during the calculation. Initially, the 
network is defined by those species with an appreciable 
abundance (> 10~ 24 mol g _1 ) plus n, p, and alphas and 
the nuclei that can be reached from any of the abundant 
species by any one of the reactions included in the net- 



work. A reaction rate is included in the network only if 
the predicted change of a molar abundance in the next 
time step, At, is larger than a threshold: 



N A p(av)Y i Y j At > 1(T 20 mol g" 



(1) 



A similar method of integration of the nuclear evolution- 
ary equations using an adaptive network has been de- 
scribed in ref. 41] . 

Our nuclear network consists of a maximum of 722 
nuclei, from free nucleons up to lnl In, linked by three 
fusion reactions: 3a, 12 C+ 12 C, and 16 0+ 16 0, electron 
and positron captures, f5~ and (3 + decays, and 12 reac- 
tions per each nucleus with Z > 6: (n, 7), (n,p), (n, a), 
(p,7)j (P' n ); (PjoOj (a,n), (a,p), (7,11), (7,p), 

and (7, a). We show the nuclear network in Table |U 
From the whole set of reactions that might be included 
in the calculations, only 3138 enter effectively into the 
reaction network equations during the integration of the 
thermodynamic trajectories in our SNIa model. 

The thermonuclear reaction rates, nuclear masses and 
partition functions are taken from the REACLIB compi- 
lation [42[ . Both theoretical and experimental thermonu- 
clear reaction rates are fitted in the JINA [55[ REACLIB 
library by an analytic function with seven parameters. 
The fits are usually better than 5% although deviations 
up to 30% are possible. The authors estimate an addi- 
tional uncertainty typically of order 30% in the original 
reaction rates. Electron screening to thermonuclear re- 
actions in the stron g, in termediate and weak regimes was 
taken into account |43l. Eij. In general, in the conditions 
achieved during thermonuclear supernova explosions the 
electron screening factors are small Weak interac- 

tion rates were taken from W9, [43 ■ 



B. Type la supernova model 

Our reference SNIa model is the one-dimensional 
delayed-detonation model DDTc in [48[ , characterized by 
its deflagration-to-detonation (DDT) transition density, 
Pddt = 2.2 x 10 7 g cm~ 3 . The supernova progenitor 
is a Chandrasekhar mass white dwarf of central density 
1.8 x 10 9 g cm -3 and uniform composition: 49.5% 12 C, 
49.5% 16 0, and 1% 22 Ne by mass. In this model the 
flame begins as a subsonic deflagration flame near the 
center of the star. As the flame propagates through the 
star, the pressure rises and the star expands. When the 
flame reaches a zone with a low enough density, Pddt, 
there is a transition to a supersonic detonation that burns 
most of the remaining fuel. Finally, the nuclear energy 
released is enough to unbind the whole star and eject its 
matter into the interstellar medium. In Fig. [1] we show 
the profiles of the most relevant physico-chemical quan- 
tities affecting the nucleosynthesis. 

This kind of SNIa model generates a layered struc- 
ture (see Fig. \2§ in which the inner several tenths of a 
solar mass achieve maximum temperatures high enough 
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TABLE I. Nuclear network 
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FIG. 1. Profiles of physico-chemical properties accross the reference model, as a function of the Lagrangian mass coordinate 
(zero at the center). Left: Peak temperature and density achieved at each mass shell during the supernova explosion (thick 
solid line). Star marks have been located every 0.1 Mq, with the center of the white dwarf at the top right end of the solid 
line, and the surface at its bottom left end. The p — T plane has been divided according to approximate locations of different 
explosive nucleosynthetic processes (dashed and dotted lines). We indicate as well the Lagrangian mass coordinate at which 
the solid line crosses the dashed and dotted lines. Right: Maximum molar fractions of neutrons, protons, and alphas achieved 
at a given mass coordinate at any time during the explosion, and final electron mole number (dot-dashed line). Note that the 
neutron molar fraction has been scaled up by a factor of 10 4 for presentation purposes. 



(Tmax ^ 5.5 x 10 9 K) to process matter into NSE, un- 
dergoing copious electron captures. When matter ex- 
pands the composition is relaxed out of NSE and consists 
mainly of iron group elements with isotopic fractions de- 
termined by the electron mole number resulting from the 
electron captures phase. As can be seen in Fig. [1] in our 
reference model the electron captures modify the progen- 
itor electron mole number only in the central ~ 0.1 M Q . 
The zone were the transition from deflagration to deto- 
nation takes place, at a mass coordinate of ~ 0.2 M Q , 
can be identified by the trough in the Y p , Y a , and Y n 
profiles. The central 0.4 Mq reach NSE, from which 



roughly 0.24 Mq experience a moderately a- rich freeze- 
out. Shortly after the detonation forms, it propagates 
fast through the white dwarf, which has no time to re- 
lax its structure before the combustion front burns most 
of the remaining fuel (this condition can be identified 
in Fig. [1] by the crowding of the star symbols between 
Tg ~ 4 and ~ 5.5). Between Lagrangian mass coordi- 
nates of ~ 0.4 and ~ 1.1 M the peak temperatures 
and densities are high enough to experience Si-burning 
and achieve quasi-statistical equilibrium (QSE) of the Fe- 
group, although this group does not achieve equilibrium 
with the Si-group. Farther out from the center, a tinier 



5 



amount of mass is subject to explosive oxygen and neon 
burning, and only a few thousandths of a solar mass ex- 
perience only explosive carbon burning. The mass of 
unburned carbon ejected by the supernova explosion is 
on the same order, in agreement with the upper limits 
deduced by [lij]. We note that all the nucleosynthetic 
processes deemed relevant in SNIa feature in our refer- 
ence model. 

For reference, we give in Table [XT] the nucleosynthe- 
sis obtained for this supernova model. The composition 
given in this and forthcoming tables corresponds to a 
time of one day after beginning of the explosion, hence 
there appear radioactive as well as stable nuclides. We 
have included in this table all nuclides whose ejected mass 
is mi > 10~ 5 M©, with the exception of 26 Al, that has 
been included because it is an interesting radionuclide. 
The ejected mass of 56 Ni is 0.675 M©, and the kinetic en- 
ergy of the ejecta is 1.16 x 10 51 erg, both values deemed 
typical for normal bright SNIa. The resulting chemical 
composition (Fig. [5]) compares well with the abundance 
stratification induced from observations of normal SNIa 
as, for instance, SN2003du (e.g., Fig. 8 in |50] , who esti- 
mated that the ejected mass of 56 Ni was 0.65 M©). Model 
DDTc also provides an excellent match to the X-ray spec- 
trum of the remnant of SN1572 (Tycho), a prototype of 
SNIa (see Fig. 7 in [HI). 
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FIG. 2. Chemical composition of the reference model as a 
function of the final velocity. The curves labelled as Fe and 
Ni include only stable isotopes. The thick curve is the mass 
fraction of 56 Ni. 

The maximum abundances of free protons, neutrons, 
and a-particles attained during the explosion are shown 
in the right panel of Fig. U as a function of the La- 
grangian mass coordinate within the exploding white 
dwarf. These profiles can be used to gain insight into the 
expected sensitivities of the nucleosynthesis with respect 



to different types of nuclear reactions, to be discussed in 
the next chapters. Neutrons are always the less abun- 
dant nucleons by ~ 5 — 6 orders of magnitude, thus we 
expect that the nucleosynthesis will not be too sensitive 
to reactions with neutrons except, perhaps, in the outer 
~ 0.1 M©. Note that neutrons are relatively abundant in 
the very center of the white dwarf, because of the lower Y e 
that results from efficient electron captures in NSE mat- 
ter at high density, but nucleosynthesis in these layers is 
not expected to be sensitive to the rate of any particular 
reaction with neutrons because the chemical composi- 
tion there is controlled by the Saha equation until mat- 
ter cools to low temperatures. Protons and a-particles 
have similar abundances within the inner ~ 1.1 M©, al- 
though their maximum molar fractions decrease steadily 
outwards within the detonated matter (M > 0.25 M©). 
Beyond ~ 1.1 M©, the maximum abundance achieved 
by protons is much lower than that of a-particles. The 
maximum temperatures attained in these layers stay be- 
low ~ 4.4 x 10 9 K, implying that the thermonuclear com- 
bustion hardly goes beyond O-burning. Thus, we expect 
that the products of O-burning will be mostly sensitive 
to reactions with a-particles. 

The above analysis can be complemented with an ex- 
amination of the molar fluxes due to different reaction 
types, e.g. (p, n), etc. In Fig. [3] we show the evolution 
of the net molar fluxes in two representative mass shells 
of our SNIa model, grouped by reaction type. The net 
molar fluxes of a given reaction type in a mass shell are 
accumulated in time according to: 



E 



(2) 



where p is density, Na is Avogadro's number, Yj is the 
molar fraction of species j, and the time integral extends 
from thermal runaway until the temperature goes below 
10 s K. The summation extends to all reactions of the 
given type, from which their inverse reactions are sub- 
tracted, e.g. in the computation of the net molar fluxes 
of the (p, n) type reactions all the (n, p) reactions are 
considered inverse reactions and their contributions are 
deducted from those of the direct, (p,n), reactions. 

In the left panel of Fig. [3] we show the net molar fluxes 
in a mass shell located at a Lagrangian mass coordinate 
of 0.5 M©. This layer was hit by the detonation wave 
~ 2.5 s after central thermal runaway, when its density 
was 1.17 x 10 7 g cm~ 3 , and heated to > 2 x 10 9 K by 
the shock front associated with the detonation. Above 
this temperature, it is the energy release by nuclear reac- 
tions which controls the evolution of temperature. The 
temperature rises very fast at the beginning due to rapid 
burning of carbon and oxygen, mainly to produce silicon 
and sulfur. About 1 ms after being shocked, a maximum 
temperature of 5.12 x 10 9 K is achieved. Later, mat- 
ter expands and cools with a longer timescale (it takes 
0.1 s to cool by 2 x 10 9 K) while most of the nuclear 
reactions are nearly in equilibrium with their inverse re- 
actions. During the heating phase, it is the 12 C + 12 C 
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TABLE II. Nucleosynthesis of the reference Type la supernova model 
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3.2-4.0 















a Range of maximum temperatures achieved in the shells in which 90% of each nuclide is produced. 



reaction which dominates the nuclear fluxes, followed by 
radiative captures of protons and a-particles once the 
temperature exceeds ~ 2.5 x 10 9 K. Compared to the 
plethora of reactions with light particles unleashed by the 
carbon fusion reaction, the contribution of the 16 O + 16 O 
reaction is quite modest until the temperature exceeds 
> 4 x 10 9 K. Above ~ 5 x 10 9 K there is a sharp in- 
crease in the cumulative molar fluxes belonging to (a, 7), 
(7,p), and (p, a) reactions, which reach similar levels. 
On the other hand, during the cooling phase there is lit- 
tle additional contribution to the net molar fluxes, and 
(p, n) and (n, a) reactions attain a level similar to that of 
16 + 16 0, while (n, 7) reactions are the ones that process 
the smallest mass. Note that, in this mass shell, the final 
cumulative molar flux due to the 16 O fusion reaction is 
about a factor four smaller than that due to the 12 C fu- 
sion reaction. Taking this mass shell as representative of 
layers that experience incomplete Si-burning, we expect 
that the products of this nucleosynthetic process will be 
most sensitive to (a, 7), (7,p), and (p, a) reactions, and 
their inverses. Note that in shells that achieve a tempera- 
ture high enough to reach NSE all the molar fluxes estab- 
lished prior to NSE are irrelevant, because NSE erases all 
memory of previous nuclear processes with the exception 
of weak interactions. 

The right panel of Fig. [3] shows the net molar fluxes in 
a mass shell located at a Lagrangian mass coordinate of 
1.2 Mq. In this case the maximum temperature achieved 
was 3.45 x 10 9 K, because the density at the time of det- 
onation impact (at t ~ 2.64 s) was only 2.4 x 10 6 g cm -3 . 
Due to the small value of the maximum temperature, the 
12 C-|- 12 C reaction dominates the molar fluxes at all times. 



Oxygen burning is incomplete, the final molar flux due 
to the 16 O fusion reaction being about 20 times smaller 
than that due to the 12 C fusion reaction. Even (n, 7) 
reactions process more matter than the 16 O fusion reac- 
tion. Taking this mass shell as representative of layers 
that do not go beyond carbon burning, we expect that 
the products of this nucleosynthetic process will be most 
sensitive to the rate of the 12 C + 12 C reaction and, to a 
lesser extent, to (p, 7) and (a, 7) reactions. 



C. Selection of the nuclear reactions 

As explained before, only 3138 nuclear reactions exceed 
the threshold of Eq. [1] and are actually included in the 
nucleosynthesis calculation. However, most of these re- 
actions contribute negligibly to the determination of the 
final chemical composition of the supernova ejecta. In 
order to determine the most relevant reactions, we define 
the total mass processed by a nuclear reaction, between 
particle k and nucleus j, in all the mass shells of the 
supernova model SDtjfc : 



m jk =J2 Ma f pWjkNtXjY,, {Aj + Ak) At 



(3) 



where M a is the mass of shell a of the supernova model, 
and Aj is the baryon number of species j. In the compu- 
tation of the integral we have not taken into account reac- 
tions above 5 x 10 9 K, because at such temperatures the 
direct and inverse reactions are in equilibrium, causing 
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Temperature (10 9 K) Temperature (10 9 K) 

FIG. 3. (Color online) Cumulative net molar fluxes of direct reactions minus inverse reactions, grouped by reaction type, 
compared to the 12 C-fusion and le O-fusion reaction fluxes, as functions of temperature. Each reaction type is identified by a 
different line type and a label, in each label an * is drawn to recall that the fluxes take into account direct and inverse reactions. 
Thin (black) curves represent those cases in which the (accumulated) molar flux of the direct reaction is larger than that of the 
inverse reaction, while thick (red) curves belong to the opposite case. Vertical lines are drawn to separate the heating phase 
from the cooling phase. Left: Evolution at Lagrangian mass coordinate of 0.5 Mq. Right: Evolution at Lagrangian mass 
coordinate of 1.2 Mq. 



the nuclear abundances to be determined by properties 
of the nuclei involved (mass, partition function) instead 
of the reaction rates. For mass shells that went through 
NSE, the computation of the integral in Eq.[3]starts when 
the temperature drops below 5 x 10 9 K, since their chemi- 
cal composition is insensitive to the nuclear history prior 
to the NSE state (with the exception of weak interac- 
tions, whose effect is not addressed in the present work). 

The reactions we have selected for careful study are 
the three fusion reactions plus those for which 9Jljk > 
10 -8 Mq. This warrants that we test all the reactions 
able to contribute significantly to the synthesis of every 
species whose yield is larger than the chosen 10 -8 Mq. 
Each time we integrate the nuclear evolutionary equa- 
tions we modify by the same factor the direct and in- 
verse reactions. Following this procedure, we find that 
the nucleosynthesis at this chosen level could be sensi- 
tive to 1096 (pairs of) reactions in addition to the above 
mentioned three fusion reactions. 

Table UlTl gives the masses processed by the three fusion 
reactions and the top ten radiative captures and trans- 
fer reactions, where the masses processed by the inverse 
reactions have been subtracted from those of the direct 
reactions. The quoted values of VJljk give a quite gen- 
erous upper limit of the impact these reactions might 
have on the resulting nucleosynthesis of the supernova, 
as the subsequent nuclear reactions destroy the products 
of earlier reactions. As we will see in the following, the 
top ten reactions listed in Table IIIII are not in fact the 



most influential reactions. 



D. Modification of the reaction rates 

As a first approach to study the sensitivity to the dif- 
ferent reaction rates, we modify them, one by one, by a 
fixed factor, either equal to f = 10 or / = 0.1, repeat- 
ing the nucleosynthesis calculation for each variation. As 
mentioned previously, each time we modify the rate of a 
reaction we modify as well by the same factor the rate 
of the inverse reaction, in order to maintain detailed bal- 
ance. 

The Gamow energies in the reactions that play a sig- 
nificant role in the nucleosynthesis of Type la super- 
novae go from a few tenths of a MeV (for instance, 
E a = 0.39 MeV for the 12 C + p <=t 13 N + 7 reaction at 
T = 10 9 K) to nearly ten MeV (e.g., E = 8.50 MeV for 
the 62 Zn + a <± 65 Ga + p reaction at T = 5 x 10 9 K). It is 
expected, both from theoretical and experimental argu- 
ments, that the uncertainties in the rates at low temper- 
atures are larger than at high temperatures 52j . Most of 
the theoretical reaction rates we have used are based on 
an statistical model of nuclei, which assumes formation of 
a compound nucleus with a high level density, a condition 
generally satisfied at high temperatures. Furthermore, 
experimental measurements of nuclear cross sections in- 
volving high-Z nuclei are generally difficult to perform at 
energies below the Coulomb barrier. Consequently, we 
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TABLE III. Masses processed by the fusion reactions and the top ten radiative captures and transfer reactions 



Reaction 


m jk (M ) 


Reaction 


3% (M Q ) 


c + c 


0.524 


M Si + a Z± "S + n 


0.70 


1Gq + 16q 


0.198 


30 P + a^ 33 S + p 


0.68 


3a 


6.7 x 10" 4 


»Si + p ^ 30p + 7 


0.67 


28 Si + aPt 32 S + 7 


0.93 


32 S + aP> 35 Cl + p 


0.65 


28 Si + p^ 29 P + 7 


0.84 


33 S + a 36 Ar + n 


0.64 


2 9p + a^ 32 S + p 


0.83 


27 Al + a^ 30 Si + p 


0.63 


28 Si + Q ^31p +p 


0.77 







use a second approach in which the reaction rates are 
modified by applying a factor that is a monotonic de- 
creasing (exponential) function of the temperature. We 
have applied the following temperature dependent factor 
to each reaction rate: 

/ (T) = l + (/ -l)ex P (-^^) ) (4) 

where f = 10 or 0.1, is the fixed factor applied in the 
first approach. Of course, we are not trying to convey 
that Eq. [4] is representative of the uncertainty of all the 
reactions studied here (see Sections IIVCI and IIVD[) , but 
it provides a convenient way to invetigate the effects of a 
temperature dependent rate error. 



amount of mass that does not experience complete car- 
bon or oxygen burning. Figure @] shows the final chemi- 
cal profiles in the outermost 0.2 M of ejecta, where the 
changes in the 12 C + 12 C and the 16 + 16 reaction rates 
are most influential. The three panels show the profiles 
belonging to our reference model and the models in which 
either the carbon or the oxygen fusion rates are increased 
by a factor of ten, both in the hydrodynamic as well as in 
the nucleosynthetic codes. As can be seen, increasing the 
12 C + 12 C rate by a factor of ten barely affects the lim- 
its of the region undergoing carbon burning, which move 
outwards ~ 0.004 M Q . On the other hand, when the 
16 O + 16 O reaction rate is enhanced by the same factor 
the limits of the oxygen burning region move outwards 
~ 0.082 M©. We conclude that the impact of the rates 
uncertainties on the energy of the supernova is negligible. 



III. SENSITIVITY TO THE RATE OF FUSION 
OF CARBON AND OF OXYGEN, AND THE 
TRIPLE-ALPHA REACTION. 

We have checked the effect of varying each fusion reac- 
tion rate by the factors given above, either taking them 
fixed or as function of the temperature. Because the fu- 
sion reactions are relevant for the nuclear energy genera- 
tion in the supernova explosion, we have recomputed the 
hydrodynamics with the modified reaction rates and give 
the results in Section UTIAI When we kept unchanged the 
thermodynamic trajectories of the reference model, but 
the reaction rates were modified in the nucleosynthetic 
code, we obtained the results shown in Section IIIIB1 



A. Rate modified in the hydrodynamic explosion 
model 

The nuclear energy release of the supernova is more 
sensitive to the rate of the 16 fusion reaction than to 
that of 12 C. The final kinetic energy of the ejecta varies 
by less than 1% when the 12 C + 12 C reaction rate is var- 
ied by a factor of 10 or 0.1, either fixed or as a function 
of the temperature given by Eq. ((4]). In contrast, the 
same relative variation in the rate of the 16 O + 16 O reac- 
tion produces a change of kinetic energy of up to ±4%. 
We ascribe this lack of sensitivity to the relatively small 



Reference model _ 




10x("0 + i 6 0) 




0.1 0.01 0.001 

Mass coordinate (M Q ) 



FIG. 4. Final chemical profile within the outer layers of the 
SNIa ejecta for three of the computed models: our reference 
model (top), the model with the 12 C + 12 C reaction increased 
by a constant factor of ten (middle), and the model with 
the le O + 16 reaction increased by a constant factor of ten 
(bottom). In this plot, the mass coordinate is zero at the 
white dwarf surface and increases inwards. 

Figures [S] to [5] and Table IIVI show the impact of the 
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FIG. 5. (Color online) Ratio of mass ejected for each nuclide 
with a modified 12 C+ 12 C reaction rate with respect to the 
mass ejected in the reference model, as a function of the mass 
fraction in the reference model (most abundant species are 
located to the right of each figure). Note that the species in- 
cluded in Table [TVl are those with an ejected mass larger than 
10 -5 Mq. Vertical lines link the results obtained for the same 
nuclide when the rate is either increased or decreased. The 
reaction rate was modified both in the hydrodynamics calcu- 
lation as well as in the nucleosynthetic code. Top: Rate mul- 
tiplied by a fixed factor, either xlO (green empty circles) or 
x 0.1 (red filled circles). Bottom: Rate multiplied by a factor 
function of temperature given by Eq. (jl}, with either fo = 10 
(green empty circles) or fo =0.1 (red filled circles). 



FIG. 6. (Color online) Ratio of mass ejected for each nuclide 
with a modified 16 0+ 16 reaction rate with respect to the 
mass ejected in the reference model, as a function of the mass 
fraction in the reference model (most abundant species are 
located to the right of each figure). Note that the species 
included in Table IIVI are those with an ejected mass larger 
than 10 -5 Mq. Vertical lines link the results obtained for the 
same nuclide when the rate is either increased or decreased. 
We modified the reaction rate both in the hydrodynamics 
calculation as well as in the nucleosynthetic code. Top: Rate 
multiplied by a fixed factor, either x 10 (green empty circles) 
or xO.l (red filled circles). Bottom: Rate multiplied by a 
factor function of temperature given by Eq. (j4|, with either 
fo = 10 (green empty circles) or fo = 0.1 (red filled circles). 



changes in the fusion rates of 12 C and of 16 O on the nu- 
cleosynthesis of the Type la supernova, when we modified 
the rates in the full supernova simulation. Figures[5]and[n] 
show the results sorted by final mass fraction of the prod- 
uct species. The mass fractions of the most abundant 
species are insensitive to the rate of fusion of 12 C. As one 
goes to smaller abundances, the scatter of the yield ratio 
is larger. Among the species with mass fraction greater 
than 0.01 there is only one nuclide that is significantly 
affected by the modification of the rate of 12 C 4- 12 C: not 
surprisingly it is 24 Mg. When the factor that modifies 
the 12 C fusion rate is a function of temperature, Eq. (Q}, 
the effect on the yields of all species is dramatically re- 
duced (bottom frame in Fig. [5]): no species experiences 
an increase larger than a factor of two in its abundance, 
and only a few species with quite small mass fractions 
(Xi < 10 -6 ) experience a reduction of more than a fac- 
tor of two in their yields when the 12 C fusion reaction 
rate is multiplied by a factor of ten. 

When the 16 0+ 16 rate is modified (Fig.|6]) the impact 
is in general smaller than when the 12 C fusion rate was 
modified. However, many of the most abundant species 
are more sensitive to the 16 O fusion rate than to the 



12 C rate because the products of C-burning (mainy 16 O, 
20 Ne, and 24 Mg) are in general less abundant than the 
products of O-burning (mainly 28 Si, 32 S, 36 Ar, and 40 Ca). 

Figure [7] presents the same results as Fig. [5] from an- 
other perspective: the impact of the modification of the 
12 C + 12 C reaction rate is shown against the element 
atomic number. The trend that can be observed in this 
figure is that increasing the 12 C fusion rate (green empty 
circles) decreases the abundances both of CNO nuclei 
and of IMEs between Phosphorus and Titanium, and in- 
creases the abundances of Magnesium, Aluminum, and 
Silicon, while elements beyond Vanadium are scarcely af- 
fected at all. If the rate of 12 C fusion is decreased (red 
filled circles) the trend is inverted, but the yields are in 
general more sensitive to a decrease in this rate than to 
an increase by the same factor. 

As can be seen in Fig. HI an increase in the 16 + 16 
reaction rate results in a small decrease in the production 
of elements up to Magnesium and an increase in elements 
from Chlorine to Chromium. The effect on the mass 
fractions is much smaller than that due to variations in 
the 12 C fusion rate. 

We give in Table ITvl the sensitivity of the yield of each 
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FIG. 7. (Color online) Same as Fig.[5]but plotted as a function 
of the atomic number of the product nucleus. Note that not 
all the isotops shown here appear in Table IIVI Top: Rate 
multiplied by a fixed factor, either x 10 (green empty circles) 
or xO.l (red filled circles). Bottom: Rate multiplied by a 
factor function of temperature given by Eq. Q, with either 
fo — 10 (green empty circles) or fo = 0.1 (red filled circles). 



one of the species included in Tablellllto the rate of fusion 
reactions. There, Di is the logarithmic derivative of the 
mass ejected of species % with respect to the enhancement 
factor of each fusion reaction, fo (note that when using 
Eq. HI fo represents the maximum enhancement factor, 
attained at low temperatures), 



Di 



d log rrii 
d login 



0.5 log 



m iA0 

mini 



(5) 



where m^io is the mass ejected of species i for fo = 10., 
and rriifi.i is the corresponding mass when f = 0.1. Ac- 
cording to this definition, a value of Di s» 0.3 means that 
the abundance of species i approximately doubles for a 
constant enhancement factor of fo = 10 in the corre- 
sponding fusion reaction rate. Similarly, a relative change 
in the abundance of a species by 12% would correspond 
to Di « 0.05, and a change by 2% would derive from 
Di « 0.01. 

Most notable is the robustness of the production of 
most Fe-group isotopes, notably of 56 Ni. When the en- 
hancement factor is computed from Eq. 21 there is no 
species with \Di\ > 0.1, neither with respect to the rate 
of 12 C + 12 C nor with respect to the 16 + 16 rate, with 
the exceptions of 12 C and 39 K, respectively. 



B. Rate modified only in the nucleosynthetic code 
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FIG. 8. (Color online) Same as Fig.[6]but plotted as a function 
of the atomic number of the product nucleus. Note that not 
all the isotops shown here appear in Table IIVI Top: Rate 
multiplied by a fixed factor, either x 10 (green empty circles) 
or xO.l (red filled circles). Bottom: Rate multiplied by a 
factor function of temperature given by Eq. (T3J), with either 
fo — 10 (green empty circles) or fo = 0.1 (red filled circles). 



Figures [9] and \W\ show the impact of the changes in 
the fusion rates of 12 C and of 16 O on the nucleosynthesis 
of Type la supernovae when the rates are modified only 
in the nucleosynthetic code. They can be compared with 
Figs. [7] and [5J respectively, to evaluate the relevance of 
incorporating the modified rates into the hydrodynamic 
code. The trends visible in these figures are qualitatively 
similar, irrespectively if the reaction rate has been mod- 
ified in the hydrodynamic calculations or not. 

In Fig. [11] we show the yield ratios belonging to modi- 
fied rate of the 3a reaction. The influence of the rate of 
this reaction focuses on a few elements: Nitrogen, Nickel, 
Copper, and Zinc (specially the isotopes of Nickel and 
Zinc produced during alpha-rich freeze-out of NSE) in- 
versely correlate with the factor of enhancement of the 
3a reaction, while Titanium and, to a lesser extent, Scan- 
dium, Manganese, and Iron (specially the isotopes pro- 
duced during explosive Si-burning) are more abundant 
when the 3a reaction is faster. These results can be ex- 
plained by the fact that for a faster rate an alpha-rich 
freeze-out occurs at lower temperature and density. As 
a result, the dotted-line in Fig. [T] (left) shifts down when 
the 3a reaction is faster, increasing the yield of species 
made in normal freeze-out at the expense of alpha-rich 
freeze-out products. 

We give in Table [V] the sensitivity of the yield of each 
one of the species included in Tablellllto the rate of fusion 
reactions, when they are modified only in the nucleosyn- 
thetic code. This table can be compared to Table IIVI 
to avaluate the importance of running a hydrodynamic 
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TABLE IV. Sensitivity of the nucleosynthesis to the rate of fusion reactions: Rate modified in the hydrodynamic and nucle- 
osynthetic codes (see also Figs. [S] to [S]o 



Nucleus 


Di ( 12 C + 12 C) 


A ( 12 C + 12 C)_ 


A ( lb O + 15 0) 


A ( O + 1£> 0)_ 




-4.8E-1 


-1.7E-1 


-3.7E-2 


-1.6E-2 


16 q 


1.9E-2 


6.0E-3 


-1.3E-1 


-5.2E-2 


20 Ne 


6.0E-2 


2.8E-2 


-5.8E-2 


-2.6E-2 


23 Na 


-3.1E-1 


-7.9E-2 


-6.4E-2 


-2.6E-2 


24 Mg 


3.7E-1 


9.2E-2 


-1.1E-1 


-4.6E-2 


25 Mg 


-2.8E-1 


-4.2E-2 


-3.0E-2 


-5.0E-3 


26 Mg 


-2.8E-1 


-5.4E-2 


-8.2E-2 


-3.7E-2 


26 Al 


-4.1E-2 


4.0E-2 


1.3E-2 


6.0E-3 


27 Al 


2.5E-1 


7.4E-2 


-6.5E-2 


-2.7E-2 


28 Si 


1.1E-2 


6.0E-3 


7.2E-2 


2.9E-2 


29 Si 


1.3E-1 


1.6E-2 


-2.5E-2 


-1.0E-3 


30 Si 


2.8E-1 


7.0E-2 


-1.6E-1 


-7.1E-2 


31p 


3.1E-2 


-1.0E-3 


-1.3E-2 


3.0E-3 


32 g 


-1.4E-2 


0. 


5.7E-2 


2.2E-2 


33 g 


-1.9E-2 


-9.0E-3 


1.7E-2 


8.0E-3 


34g 


-6.7E-2 


-5.0E-2 


-7.2E-2 


-3.0E-2 


35 CI 


-2.5E-1 


-6.7E-2 


1.2E-1 


6.0E-2 


36 Ar 


-1.9E-2 


-2.0E-3 


1.6E-2 


5.0E-3 


37 Ar 


-1.5E-1 


-3.0E-2 


1.8E-1 


7.8E-2 


38 Ar 


-1.9E-1 


-4.2E-2 


1.2E-1 


4.8E-2 


39 K 


-2.2E-1 


-3.6E-2 


2.3E-1 


1.0E-1 


40 Ca 


-2.3E-2 


-6.0E-3 


-2.1E-2 


-1.1E-2 


44 Ti 


-3.1E-2 


-9.0E-3 


-5.6E-2 


-2.8E-2 


48 y 


-2.4E-2 


-1.0E-2 


-4.1E-2 


-1.6E-2 


49 y 


-1.9E-2 


-5.0E-3 


-2.2E-2 


-1.1E-2 


50 Cr 


-4.0E-3 


6.0E-3 


8.0E-2 


3.2E-2 


51 Cr 


-1.2E-2 


-3.0E-3 


-9.0E-3 


-4.0E-3 


52 Mn 


-2.0E-2 


-1.2E-2 


-4.0E-2 


-1.5E-2 


53 Mn 


-1.4E-2 


-7.0E-3 


-2.8E-2 


-1.2E-2 


54 Fe 


0. 


3.0E-3 


1.3E-2 


5.0E-3 


55 Fe 


0. 


0. 


0. 


0. 


56 Fc 
r e 


u. 


9 OF 3 


^ OF 3 


O.U-C/-0 


57 Co 


-5.0E-3 


-4.0E-3 


-1.3E-2 


-5.0E-3 


56 Ni 


-9.0E-3 


-9.0E-3 


-2.8E-2 


-1.2E-2 


58 Ni 


2.0E-3 


0. 


0. 


0. 


59 Ni 


1.0E-3 


0. 


0. 


0. 


60 Ni 


2.0E-3 


-1.0E-3 


1.0E-3 


0. 


61 Ni 


4.0E-3 


2.0E-3 


2.0E-3 


-2.0E-3 


62 Ni 


6.0E-3 


2.0E-3 


0. 


-2.0E-3 



a Values of Di less than 1.0E-3 have been put to 0. 

b Enhancement factor function of temperature according to Eq. [4] 



code with the reaction rates modified or take the thermo- 
dynamic profiles of a reference model and modifying the 
rates only in a post-processing code. One finds that the 
sensitivities shown in both tables are qualitatively simi- 
lar. Although the precise values of Di for given species 
are not equal, the rating of the species that are most 
sensitive to any fusion reaction rate is the same in both 
tables. Given the pre-eminence of the fusion reaction 
rates with respect to the release of nuclear energy, we 
conclude that for SNIa this kind of study can be safely 
carried out with a post-processing code, using a set of 
thermodynamic trajectories obtained with a supernova 
hydrodynamics code where the reaction rates remain un- 



changed. 



IV. SENSITIVITY TO THE RATE OF 
RADIATIVE CAPTURES AND TRANSFER 
REACTIONS 

We will discuss in this section the sensitivity of the nu- 
cleosynthesis to changes in the rate of radiative captures 
and transfer reactions. We measure the sensitivity in a 
similar way as with respect to the rate of the fusion re- 
actions, by defining Di as in Eq. [S] The meaning of Di 
is now the logarithmic derivative of the mass ejected of 
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TABLE V. Sensitivity of the nucleosynthesis to the rate of fusion reactions: Rate modified only in the nucleosynthetic code 
(see also Figs. [gito[TT|PI. 



Nucleus 


Di ( 12 C + 12 C) 


Di ( 12 C + 12 C)^ 


A ( O + ls O) 


Di ( ie O + 16 0)^ 


Di (3a) 


Di (3q^ 




-3.8E-1 


-1.9E-1 


0. 


0. 


0. 


0. 


16q 


-1.3E-2 


-6.0E-3 


-7.2E-2 


-3.1E-2 


-1.0E-3 


0. 


20 Ne 


1.0E-1 


5.2E-2 


0. 


0. 


0. 


0. 


23 Na 


-1.7E-1 


-6.5E-2 


-1.3E-2 


-5.0E-3 


-1.0E-3 


0. 


24 Mg 


3.9E-1 


1.1E-1 


-6.8E-2 


-3.2E-2 


-1.0E-3 


0. 


25 Mg 


-6.4E-2 


-2.0E-2 


3.0E-2 


1.7E-2 


0. 


0. 


26 Mg 


-1.1E-1 


-4.3E-2 


-4.3E-2 


-2.1E-2 


-1.0E-3 


0. 


26 Al 


2.4E-2 


5.0E-3 


7.5E-2 


4.3E-2 


0. 


0. 


27 A1 


2.5E-1 


7.4E-2 


-3.0E-2 


-1.4E-2 


-1.0E-3 


0. 


28 Si 


5.0E-3 


0. 


7.4E-2 


3.4E-2 


-2.4E-2 


-7.0E-3 


29 Si 


1.4E-1 


2.7E-2 


7.9E-2 


4.3E-2 


-2.0E-3 


-1.0E-3 


30 Si 


2.8E-1 


8.5E-2 


-1.2E-1 


-6.4E-2 


-2.0E-3 


-1.0E-3 


31p 


-1.3E-2 


-1.6E-2 


6.0E-2 


3.2E-2 


-2.0E-3 


0. 


32 g 


-2.3E-2 


-7.0E-3 


3.6E-2 


1.5E-2 


-1.7E-2 


-5.0E-3 


33 g 


-9.6E-2 


-4.1E-2 


9.0E-2 


4.3E-2 


-1.0E-3 


0. 


34 g 


-9.6E-2 


-7.0E-2 


1.1E-2 


7.0E-3 


-2.0E-3 


-1.0E-3 


35 CI 


-3.2E-1 


-9.5E-2 


1.7E-1 


8.7E-2 


1.0E-3 


3.0E-3 


36 Ar 


-2.3E-2 


-7.0E-3 


-2.0E-2 


-9.0E-3 


-7.0E-3 


-2.0E-3 


37 Ar 


-2.1E-1 


-5.5E-2 


1.7E-1 


7.0E-2 


2.0E-3 


0. 


38 Ar 


-2.4E-1 


-7.4E-2 


1.4E-1 


6.9E-2 


-2.0E-3 


0. 


39 K 


-2.7E-1 


-6.0E-2 


1.7E-1 


7.6E-2 


1.7E-2 


1.0E-2 


40 Ca 


-2.6E-2 


-9.0E-3 


-7.2E-2 


-3.3E-2 


6.0E-3 


2.0E-3 


44 Ti 


-2.4E-2 


-9.0E-3 


-9.2E-2 


-4.2E-2 


1.4E-1 


9.6E-2 


48 y 


-7.0E-3 


-2.0E-3 


-6.8E-2 


-3.0E-2 


2.9E-2 


1.5E-2 


49 y 


-8.0E-3 


-1.0E-3 


-4.1E-2 


-1.7E-2 


5.0E-3 


0. 


50 Cr 


-6.0E-3 


3.0E-3 


5.6E-2 


2.4E-2 


1.4E-2 


5.0E-3 


51 Cr 


0. 


3.0E-3 


-1.8E-2 


-7.0E-3 


1.7E-2 


5.0E-3 


52 Mn 


2.0E-3 


-1.0E-3 


-7.3E-2 


-3.1E-2 


3.5E-2 


1.3E-2 


53 Mn 


2.0E-3 


3.0E-3 


-4.4E-2 


-1.8E-2 


3.2E-2 


1.0E-2 


54 Fe 


2.0E-3 


3.0E-3 


3.0E-3 


2.0E-3 


3.5E-2 


1.3E-2 


55 Fe 


2.0E-3 


3.0E-3 


-2.6E-2 


-1.2E-2 


4.8E-2 


1.7E-2 


56 Fe 


n 


n 
u. 




n 


n 

U. 


n 
u . 


57 Co 


2.0E-3 


0. 


-6.0E-3 


-2.0E-3 


-1.1E-2 


-4.0E-3 


56 Ni 


3.0E-3 


0. 


-1.6E-2 


-8.0E-3 


1.2E-2 


5.0E-3 


58 Ni 


2.0E-3 


0. 


-1.0E-3 


0. 


-4.7E-2 


-1.7E-2 


59 Ni 


0. 


0. 


0. 


0. 


-5.6E-2 


-2.2E-2 


60 Ni 


2.0E-3 


0. 


0. 


0. 


-1.1E-1 


-4.5E-2 


61 Ni 


2.0E-3 


0. 


0. 


0. 


-1.4E-1 


-5.8E-2 


62 Ni 


2.0E-3 


0. 


0. 


0. 


-1.6E-1 


-6.4E-2 



a Values of Di less than 1.0E-3 have been put to 0. 

b Enhancement factor function of temperature according to Eq. [4] 



species i with respect to the enhancement factor, /o, of a 
reaction between particle k and nucleus j, while m^io is 
the mass ejected of species i when an enhancement fac- 
tor /o = 10. (either fixed or function of temperature) is 
applied to the rate of the reactions j + k — > I + m and 
I + m — > j + k, and mi o i is the corresponding mass when 
/o = 0.1. 

We start by analyzing the results obtained with a fixed 
enhancement factor (our first approach). In Section llVBI 
we present the results obtained when we compute the 
enhancement factor with a decreasing uncertainty (see 
Eq. 21 our second approach) . Then, we select the reac- 
tions to which the nucleosynthesis is most sensitive and 



analyze in Section IIV CI the results achieved by adopt- 
ing different prescriptions for their reaction rates, cho- 
sen among the most recent literature. Finally, in Sec- 
tion IIVDI we analyze the temperature ranges in which 
a modification of a given reaction rate affects most the 
chemical composition of the supernova ejecta. 



A. Fixed rate enhancement factor 

Tables I VII to IXII give, for each reaction pair j + k £± 
I + m that has a significant impact on the nucleosynthe- 
sis, the nuclei i for which |Dj| > 0.3 (more than twofold 
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FIG. 9. (Color online) Ratio of mass ejected for each nu- 
clide with a modified 12 C+ 12 C reaction rate with respect to 
the mass ejected in the reference model, as a function of the 
atomic number of the product nucleus. We modified the re- 
action rate only in the nucleosynthetic code. Top: Rate mul- 
tiplied by a fixed factor, either xlO (green empty circles) or 
x 0.1 (red filled circles). Bottom: Rate multiplied by a factor 
function of temperature given by Eq. @, with either fo = 10 
(green empty circles) or fo =0.1 (red filled circles). 



FIG. 10. (Color online) Ratio of mass ejected for each nu- 
clide with a modified 16 0+ 16 reaction rate with respect to 
the mass ejected in the reference model, as a function of the 
atomic number of the product nucleus. We modified the re- 
action rate only in the nucleosynthetic code. Top: Rate mul- 
tiplied by a fixed factor, either x 10 (green empty circles) or 
xO.l (red filled circles) . Bottom: Rate multiplied by a factor 
function of temperature given by Eq. (T3J), with either fo — 10 
(green empty circles) or fo = 0.1 (red filled circles). 



increase or decrease in the yield when the rate in en- 
hanced or decreased by a factor of ten), and those for 
which 0.3 > \Di\ > 0.05 (relative increase or decrease of 
the yield between 12% and a factor of two). Although 
we only list in the tables the direct reactions, the inverse 
reactions contributed as well to the changes in the nu- 
cleosynthetic yield. There is only one reaction pair for 
which \Di\ > 1, it is the 30 Si + p £± 31 P + 7 reaction 
and the species whose abundance is mostly affected is 
35 S. Each table shows the reactions belonging to a given 
type, e.g. (n, 7), sorted according to the total mass they 
processed in our reference model, fSljk (see Table ITTT|) . 

The species most sensitive to changes in the rates 
of (n, 7 ) reactions (Table Ell) are 17 0, 26 Al, 21 Ne, and 
35 S. The yields of all these nuclides are small, of order 
10~ 7 Mq. Apart from the neutron captures on iron iso- 
topes, all the reactions listed in Table I VII involve IMEs 
or CNO elements as the parent nuclides. Among the 
species with \Di\ > 0.05, the most abundant are 29 Si 
and 33 S, both with yields on the order of a few times 
10 -4 Mq, suggesting that the temperature range where 
(n, 7) reaction rates most affect the final nucleosynthesis 
is approximately 2 x 10 9 <T<4x 10 9 K, in agreement 
with our analysis in Section III Bl 

Radiative captures of protons are by far the group of 
reactions whose rate most strongly determine the final 
abundances of the supernova explosion, as can be de- 



duced from Table \VU\ The reaction with the largest \Di\ 
of the whole network is 30 Si + p 31 P + 7, for which as 
many as 20 product species have \Di\ > 0.05. The species 
most affected by changes in the rates of the proton cap- 
ture group of reactions are 14 N, 26 Mg, 26 Al, 27 A1, 32 P, 
35 S, and 43 Ca. Among these, 27 A1 is the species with the 
largest yield, m s = 4.6 x 10~ 4 M . Within the species 
with 0.3 > \Di\ > 0.05 there are important products of 
the supernova explosion such as 24 Mg, 25 Mg, 29 Si, 30 Si, 
31 P, 33 S, 34 S, 35 C1, 38 Ar, 50 Cr, 51 Cr, and 55 Fe. The par- 
ent nuclei involved in these reactions cover a wide range 
from 14 C to 64 Ga. 

The rates of (p, n) reactions do not influence signif- 
icantly the nucleosynthesis. The most affected species 
are 32 P and 37 C1, both with final yields on the order of 
10 -7 Mq. Among the species with |D,-| > 0.05, the most 
abundant are 29 Si, with = 4.5 x 10 -4 Mq, 25 Mg and 
49 V, both with yields on the order of 10~ 5 M . 

There are several reactions of radiative capture of a 
particles that bear a non-negligible influence on the syn- 
thesis of large numbers of species (Table IIX|) . The most 
notable is the reaction 20 Ne + a 24 Mg + 7, for which 
there are 13 species with \Di\ > 0.3 and 20 species with 
0.3 > \Di\ > 0.05, among them species with large abun- 
dances such as 20 Ne, 24 Mg, 28 Si, 36 Ar, 40 Ca, or 52 Mn. 
Variations in the rate of (a, 7) reactions on 24 Mg and 
12 C also influence the yields of large numbers of species, 
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TABLE VI. Sensitivity of the nucleosynthesis to the rate of (n, 7) reactions with the parent nuclide given in the first columnQ 



Parent nuclide 


m, k (M ) 


Nuclei with |A| > 0.3 


Nuclei with 0.3 > |A| > 0.05 


28 Si 


0.30 




2i Ne, 2b Mg, 2b Al, 33 S 


55 Fe 


0.20 




54 ' 55 Mn, 57 > 58 Fe, 58 ' 59 Co 


32 g 


0.17 


35 s 


25 Mg, 26 Al, 29 Si, 32 ' 33 P, 33 S 


36 Ar 


0.093 




37 CI, 37 Ar 


44 Ti 


0.076 




45 Sc 


24 Mg 


0.045 


26 Al, 35 S 


17 0, 21 Ne, 25 Mg, 32 P, 33 S 


25 Mg 


0.010 


21 Ne 


17 , 25 Mg, 26 Al 


5e Fe 


0.0093 




57 Fe 


16 Q 


0.0092 


17 


25 Mg 


46 Ti 


0068 




46 Ti 


29 Si 


0.0061 




29 Si 


20 Ne 


0.0049 


21 Ne 




33g 


0.0033 




33p 


35 C1 


0.0029 




37 CI 


12 C 


0.0024 




21 - 22 Ne, 25 Mg, 45 Sc 


31p 


9.1 x 10~ 4 




32 p 



a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table ITTTT) and with any 
max(|D,|) > 0.05. 
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FIG. 11. (Color online) Ratio of mass ejected for each nuclide 
when the 3a reaction rate is modified to the mass ejected in 
the reference model, as a function of the atomic number of the 
product nucleus. We modified the reaction rate only in the 
nucleosynthetic code. Top: Rate multiplied by a fixed factor, 
either xlO (green empty circles) or xO.l (red filled circles). 
Bottom: Rate multiplied by a factor function of temperature 
given by Eq. Q, with either fo — 10 (green empty circles) or 
fo — 0.1 (red filled circles). 



although not as much as does the reaction 20 Ne(a,7), 
while captures on 28 Si and 32 S have a much more limited 
reach. Note that, because we always modify coherently 
direct and inverse reaction rates to maintain detailed 



balance, the sensitivity to the reactions 20 Ne (a, j) 24 Mg 
and 24 Mg (a, 7) Si also pick up the effect of varia- 
tions on 24 Mg (7, a) 20 Ne and 28 Si(7, a) 2 Mg, respec- 
tively, which are reactions relevant for silicon burning. 
Most of the parent nuclei listed in the table belongs to 
the IMEs group. An interesting exception is the reac- 
tion 58 Ni + a ^ 62 Zn + 7, which plays a relevant role 
in the alpha-rich freeze-out of incinerated matter that 
leaves NSE at densities below ~ 10 8 g cm -3 , because of 
the adiabatic expansion of the ejecta. 

The species most affected by variations on the rate of 
(a, n) reactions, as well as the parent species listed in 
Table [Xj belong to the IMEs group. Among the species 
with \Di\ > 0.3, the most abundant is 26 Mg, whose yield 
is 2.7 x 10~ 5 M . 

In Table IXI1 we give the sensitivities to the rates of 
(ct,p) reactions. Abundant species most affected by vari- 
ations of this type of reactions are 23 Na, 26 Mg, and 27 A1, 
the last with a yield of 4.6 x 10~ 4 M Q . As with (p, 7) 
reactions, the parent species cover a wide range of baryon 
numbers, from 13 N to 62 Zn. The most influential reac- 
tions are 13 N + a v± 16 + p, 20 Ne + a <=* 23 Na + p, 
23 Na + a ^ 26 Mg + p, 24 Mg + a V± 27 Al + p, and 
27 A1 + a <=^ 30 Si + p. There are present as well several 
reactions relevant for the alpha-rich freeze-out of NSE, 
such as 56 Ni + a 59 Cu + p and 58 Ni + a 61 Cu + p. 

Several reactions appear in the tables that are respon- 
sible for bridging the gap between the quasi-statistical 
equilibrium groups (the QSE groups of silicon and iron) 
in silicon burning [5J|, such as 45 Sc + p ^ 46 Ti + 7, 
42 Ca+a <=* 46 Ti+7, and 44 Ti+a <=► 47 V+p. It is remark- 
able that the set of abundances significantly affected by 
the modification of the rates of these reactions is limited 
to species that fall into the QSE groups gap, but there 
appear very few species belonging to the QSE groups. 
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TABLE VII. Sensitivity of the nucleosynthesis to the rate of (p, 7) reactions with parent nuclide given in the first column. 



Parent nuclide 


m Jk (m s ) 


Nuclei with A > 1. 


Nuclei with 1 > |A| > 0.3 


Nuclei with 0.3 > |A| > 0.05 


29 Si 


0.67 








21 Ne, 2b Mg, 2b Al, 3b S, 43 Ca, 4 'Ti 


57 Co 


0.63 








54,55 Mn] 58p e 


58 Ni 


0.61 








63 Cu 


34g 


0.61 








35g 


53 Mn 


0.58 








52 Cr 


35 C1 


0.56 








35 S, 37 CI, 37 Ar 


30 Si 


0.56 




35 g 


32 p 


21 Ne, 23 Na, 24 " 26 Mg, 26 Al, 29 < 30 Si, 

31,33p 33,34g 35,37qj 38 Ar 42 ' 43 Ca 
47 T j ' 


27 Al 


0.53 






27 Al, 35 S, 43 Ca 


23 Na, 24 " 26 Mg, 26 Al, 32 P, 34 S, 35 C1, 










47 Ti 


55 Co 


0.47 








50 Cr, 55 Fe, 56 Co 


39 K 


0.39 








39 K, 41 > 43 Ca 


59 Cu 


0.29 








59,60 Nii 63 CU) 64,65 Zn 


56 Co 


0.28 








56 Co 


48y 


0.27 








46,47rpj 


33g 


0.27 








43 Ca, 47 Ti 


30 p 


0.25 








43 Ca 


51 Mn 


0.22 








50,51 Cr 


47y 


0.14 








46,47rpj 


47 Ti 


0.13 








46,47rpj 


26 Al 


0.12 






26 Al 




54 Mn 


0.087 








54 Mn 


44 Sc 


0.082 








47 T j 


59 Co 


0.079 








58 Fe, 59 Co 


55 Mn 


0.065 








55 Mn 


26 Mg 


0.061 






26 Mg 


23 Na, 25 Mg, 26 Al, 29 Si, 35 S , 43 Ca, 










47 Ti 


36 CI 


0.026 








37 C1 


25 Mg 


0.023 






26 Al 


21 Ne, 25 Mg, 35 S 


43 Sc 


0.021 








43 Ca 


23 Na 


0.014 








45 Sc 


58 Fe 


0.011 








58 Fe 


45 Sc 


0.0081 








42 Ca, 45 Sc, 46 Ti 


62 Zn 


0.0068 








63 Cu 


62 Cu 


0.0060 








63 Cu, 65 Zn 


42 Sc 


0.0042 








43 Ca 


37 C1 


0.0017 








37 C1 


21 Ne 


5.4 x 10' 


-4 






21 Ne 


17 F 


3.9 x 10" 


-4 




14 N 




14 c 


1.9 x 10" 


-4 






21 Ne 


64 Ga 


6.3 x 10" 


-5 






63 Cu, 65 Zn 


63 Ga 


5.0 x 10" 


-5 






64 Zn 



a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table ITTT1) and with any 
max (|A|) > 0.05. 



When one of these reactions is modified, the increase or 
decrease it produces in the flux between QSE groups is 
offset against a slight adjustment in the abundances of 
other species located within the gap, therefore leaving 
the abundances of the species in QSE nearly unchanged. 
For instance, modifying the rate of 42 Ca + a 46 Ti + 7 
by a factor of xlO, the global flux carried by all the re- 
actions that bridge the QSE gap changes by less than 
5%, in spite of an increase in the flux carried by the 
42 C a + a <=£ 46 Ti + 7 reaction by a factor of x7 (to- 



gether with a decrease of the final yield of Ca by 
a factor of x0.7 and an increase of the final yield of 
46 Ti by a factor of xl.8). The larger flux carried by 
42 C a + a ^ 46 Ti + 7 is offset by a decrease in the fluxes 
due to 44 Ti + a <^ 47 V + p, and other reactions within 
the gap. 

We have plotted in Fig. [12] the reactions with the 
largest max(|Dj|) in Tables IVT1 to IXll It is remarkable 
that no reaction appears involving the main products 
of SNIa, i.e. elements from the Fe-group nuclei, except 
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TABLE VIII. Sensitivity of the nucleosynthesis to the rate of (p, n) reactions with parent nuclide given in the first column. 



Parent nuclide 


OT jfe (M ) 


Nuclei with |A| > 0.3 


Nuclei with 0.3 > |A| > 0.05 


M Mn 


0.61 




b3 Mn 


"Co 


0.57 




57 Co, 61 Ni, 65 Zn 


30 Si 


0.54 




25 Mg, 26 Al, 29 Si, 32 P, 35 S 


32p 


0.44 


32 p 




56 Fe 


0.31 




55 Mn, 57 - 58 Fe 


58 Co 


0.27 




58 Co 


49y 


0.26 




49y 


47 Ti 


0.12 




46 Ti 


27 A1 


0.12 




26 Al 


54 Mn 


0.11 




54 Mn 


44 Sc 


0.072 




47 Ti 


55 Mn 


0.068 




55 Mn 


45 Sc 


0.064 




45 Sc 


37 C1 


0.054 


37 CI 




62 Cu 


0.015 




63 Cu, 65 Zn 


26 Mg 


0.014 




26 Al 


18 Ne 


7.1 10~ 4 




14 N 


35 S 


3.4 10 -5 




35 g 



a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table ITTTT) and with any 
max (| |) > 0.05. 



TABLE IX. Sensitivity of the nucleosynthesis to the rate of (a, 7) reactions with parent nuclide given in the first column. 



Parent nuclide 


m jk (m ) 


Nuclei with |A| > 0.3 


Nuclei with 0.3 > |A| > 0.05 


28 Si 


0.93 




SUgj 31,33p 33,34g 35q] 38 ^ 


32 g 


0.39 




37 ci 


20 Ne 


0.33 


24,25 Mg] 26,27 ^ 30^ 35,37^ 


14 N, 20 - 21 Ne, 23 Na, 26 Mg, 28 > 29 Si, 32 < 33 P, 33 ' 34 S, 36 " 38 Ar, 






39 K, 41 " 43 Ca, 46 ' 47 Ti 


40 Ca, 44 Ti, 48 - 49 V, 52 ' 53 Mn 


16 


0.30 


14 N, 20 > 21 Ne, 23 Na, 24 Mg, 27 Al, 29 Si, 32 P, 33 ' 35 S, 45 Sc 


40 Ca 


0.20 




44 Ti 


24 Mg 


0.19 


24 Mg, 35 S 


23 Na, 25 Mg, 26 ' 27 A1, 30 Si, 31 P, 35 C1, 38 Ar, 39 K, 41 ' 42 Ca, 
45 Sc, 46 - 47 Ti 


58 Ni 


0.15 




62 Ni, 63 Cu, 64 " 66 Zn 


57 Ni 


0.090 




61 Ni 


12 C 


0.074 


45 Sc 


14 N, 28 - 30 Si, 32 P, 37 Ar, 39 K, 40 - 42 Ca, 44 - 46 Ti, 48 V, 
52 Mn 


29 Si 


0.065 




33g 


33 g 


0.062 




37 C1 37 Ar 


30 Si 


0.047 


35g 


30 Si,' 31 - 33 P, 34 S, 35 C1, 37 > 38 Ar, 39 K, 42 - 43 Ca 


34g 


0.029 




38 Ar, 39 K 


41 Ca 


0.024 




43 Ca, 47 Ti 


42 Ca 


0.011 




38 Ar, 39 K, 41 ' 42 Ca, 46 ' 47 Ti 


14 N 


3.1 10" 4 




21 Ne 


62 Zn 


1.2 10" 4 




66 Zn 


17 o 


1.0 10" 4 




21 Ne 



a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table [TTTTl and with any 
max (ID;!) > 0.05. 



47 Cr + a 50 Mn + p. This reaction might play a role 
in the freeze-out from incomplete silicon burning [54j , as 
one of the last links between the two main QSE groups. 
The rest of the reactions in the plot sketch a connected- 
network from 12 C up to 37 Ar. The pattern displayed in 
the figure suggests that we can talk not only of the re- 
actions whose rates are most influential in the supernova 



yields, but also of a path in the nuclear chart that is 
most influential. The main stream involves (a, 7) reac- 
tions from 12 C to 28 Si, although the step from 16 O to 
20 Ne is not covered by the 16 O + a ^ 20 Ne + 7 reac- 
tion pair but by the combination of 16 17 + 7 
and 17 + a <=± 20 Ne + n. Above Mg, there appear 
many branches due to a number of (a, p) and (a, n) re- 
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TABLE X. Sensitivity of the nucleosynthesis to the rate of (a, n) reactions with the parent nuclide given in the first column. 



Parent nuclide 


OT jfe (M ) 


Nuclei with |Aj > 0.3 


Nuclei with 0.3 > |A| > 0.05 


29 Si 


0.70 


3b g 


aj -Ne, 25 Mg, " b Al, 29 Si, 34 S, 3b Cl, 43 Ca, 4Y Ti 


33g 


0.64 




37 C1, 43 Ca 


27 A1 


0.54 


43 Ca, 47 Ti 


25 Mg, 26 - 27 Al, 29 Si, 32 < 33 P, 33 ' 35 S, 37 C1 


30 Si 


0.43 


33 P, 43 Ca 


29 Si, 32 P, 33 < 35 S, 37 C1, 47 Ti 


25 Mg 


0.28 


26 Al 


17 0, 21 Ne, 25 Mg, 32 P, 35 S 


34g 


0.10 


37 C1 


37 Ar 


26 Mg 


0.075 


26 Mg, 35 S 


21 Ne, 25 Mg, 26 - 27 Al, 29 Si, 32 P, 43 Ca, 47 Ti 


46 Ti 


0.058 




46 Ti 


41 Ca 


0.055 




43 Ca, 47 Ti 


38 Ar 


0.042 




41 Ca 


21 Ne 


0.014 


21 Ne 




22 Ne 


0.013 




17 0, 22 Ne 


17 


0.011 


17 o 


21 Ne 


14 C 


7.5 x 10" 4 




17 



a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table ITTTll and with any 
max (ID, |) > 0.05. 



actions and their inverses, which shift the stream to the 
side of moderately neutron-rich nuclei. The path ends 
in a loop involving the reactions 34 S + a 37 Ar + n, 

37 Ar + n ^ 37 C1 + pj and 37 C1 + p ^ 34g + a 

The right panel of Fig. [12] shows the path defined by 
the reactions that carry the largest mass flow during the 
SNIa explosion, computed through Eq. [3] together with 
the path of most influential reactions. A couple of points 
have to be retained to make a meaningful comparison be- 
tween both paths. First, as explained in Section [II C[ in 
shells achieving NSE, the mass flow prior to equilibrium 
has not been accounted for. Hence, Fig. Q~S] does not re- 
flect the nuclear flow from carbon and oxygen up to the 
Fe-group in incinerated shells, giving the impression that 
the Si-group and the Fe-group are disconnected. Second, 
our consideration of a nuclear reaction as highly influen- 
tial is based on the relative variation of the yield of any 
species that has max(|D,*|) > 0.3. Thus, it is possible 
for a nuclear reaction off the maximum mass flow path 
to affect significatively the abundance of trace species. 
The maximum mass flow path follows the Z = N line 
from carbon to the Si-group, but deviates to the neutron 
rich side within the Fe-group, that accumulates the neu- 
tron excess due to the initial metallicity and the electron 
captures close to the center of the star (see Fig.Q]). Con- 
sequently, and due to the small abundance of neutrons, 
the connection between both groups is provided mainly 
by (a, p) reactions. It is evident from Fig. [12] that the 
most influential path occupies a region slightly more neu- 
tron rich as compared to the maximum mass flow path, 
especially below A ~ 28. The species most affected by 
reactions that lay off the maximum mass flow path are 
trace species, e.g. 17 0, 21 Ne, 25 Mg, 37 C1, and 47 Ti. Reac- 
tions within the Si-group form a dense network in which 
the rate of a particular reaction loses relevance, and the 
same applies to reactions within the Fe-group. 

Tables IXIll and IXI 1 1 1 gi ve . for the most important prod- 



uct nuclei, i, the reactions j + k ^ I + m that have the 
largest impact on its yield, |D,-|. We have included in 
these tables only the species with mass fraction greater 
than 10~ 5 , or that are interesting radioactive isotopes, 
and with max(|Dj|) > 0.01. We show as well the pro- 
duction factor of these species (in our reference model), 
to help in evaluating the relevance of the modifications 
to the yield of each nuclide. For radioactive isotopes, we 
have calculated the production factor taking as a refer- 
ence the solar abundance of the end product of the disin- 
tegration chain. For each nucleus, we show a maximum 
of ten reactions, sorted by \Di\. 

Among the species with the largest production factors, 
the yields of 28 Si and 32 S are hardly affected by any rate 
(Table [XII]), the maximum \D t \ being 0.09 and 0.04 re- 
spectively (both due to the reaction 13 N + a 16 + p), 
implying relative variations on their yields of 23% and 
10% when the rates change by a factor of ten. The 
same applies to 36 Ar, whose maximum \Di\ is 0.09 (Ta- 
ble [Xiri|, while the yield of 40 Ca is slightly more depen- 
dent on the rates of the radiative a captures on 20 Ne 
and 12 C and on the (a, p) reaction on 13 N, with \Di\ up 
to 0.15 (variation of up to 40% of the yield for a rate 
change by a factor of ten). Isotopes belonging to the 
Fe-group with production factor larger than 100 have a 
similar level of sensitivity to the variation of the reaction 
rates with maximum \Dj\ slightly above ~ 0.1, with the 
exception of 56 Ni, 58 Ni, 54 Fe, and 56 Fe (the last two do 
not appear in the table), whose yields are quite robust. 

The list of reactions that are most influential on 44 Ti 
synthesis has some points in common with that found 
by [19| in the context of core-collapse supernovae as, for 
instance, 44 Ti + a T=t 47 V + p, 40 Ca + a <± 44 Ti + 7, 
12 C + a <=^ 16 + 7, and the 3a reaction. However, 
in SNIa there are no reactions involving nuclei heavier 
than Ti that affect significantly the yield of 44 Ti, at vari- 
ance with what was found by [l9| for core-collapse su- 
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TABLE XI. Sensitivity of the nucleosynthesis to the rate of (a, p) reactions with parent nuclide given in the first column. 



Parent nuclide 


OT jfc (M ) 


Nuclei with |Aj > 0.3 


Nuclei with 0.3 > |A| > 0.05 


28 Si 


0.77 




2b Mg, ^-^p, 30 S, 3 'C1 


30p 


0.68 




43 Ca 


32g 


0.65 




32,33p 35g 35,37qj 


27 Al 


0.63 


27 A1, 35 S 


24 " 26 Mg, 26 Al, 29 ' 30 Si, 32 ' 33 P, 33 S, 37 C1, 43 Ca, 47 Ti 


31p 


0.61 


30 Si 31_33 P 34,35g 35q 


24 Mg 


0.58 


35 S, 43 Ca 


24-26 Mg] 26,27 A1> 30^ 32^ 34^ 35,37^ 47 Ti 


56 Ni 


0.58 




43 Ca, 47 Ti, 59 - 60 Ni, 63 Cu, 64 ' 65 Zn 


33g 


0.55 




43 Ca 


39 K 


0.48 




43 Ca, 47 Ti 


40 Ca 


0.44 




43 Ca 


29 Si 


0.44 




32p 


13 N 


0.39 




14 N 28 gi 37,38 A 40,43 c 45g 44^ 48,49 V 5 c 52,53 Mn 

1 1 1 1 1 1 11 


35 CI 


0.38 




38 Ar, 39 K, 41 ' 42 Ca, 47 Ti 


20 Ne 


0.37 


23 Na 


17 0, 21 Ne, 26 Mg, 35 S, 43 Ca 


25 Mg 


0.24 




26 Al 


58 Ni 


0.23 




62 Ni, 63 Cu, 64 < 66 Zn 


44 Ti 


0.18 




14 N, 45 Sc, 44 ' 47 Ti, 48 - 49 V 


57 Ni 


0.17 




61 Ni, 65 Zn 


48 Cr 


0.13 




48y 49y 


45 Ti 


0.12 




45 Sc 


23 Na 


0.12 


26 Mg, 43 Ca 


14 N, 21 Ne, 23 Na, 29 Si, 32 P, 33 S, 37 CI, 40 Ca, 45 Sc, 44 - 47 Ti 


41 Ca 


0.068 




43 Ca, 47 Ti 


46 Ti 


0.068 




46 T j 


35 Ar 


0.061 




43 Ca 


34g 


0.054 


37 C1 




39 Ca 


0.049 




47 Ti 


30 Si 


0.032 




33p 


42 Ca 


0.020 




46 Ti 


60 Zn 


0.0077 




63 Cu, 64 Zn 


62 Zn 


0.0054 




66 Zn 


42 Sc 


0.0024 




43 Ca 


61 Zn 


0.0020 




65 Zn 


17 Ne 


7.6 10" 4 


14 N 




43 Ti 


4.1 10" 4 




47 Ti 


18 Ne 


3.4 10" 4 




14 N 


21 Ne 


1.6 10" 4 




21 Ne 


47 Cr 


8.9 10" 5 


47 Ti 





a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table ITTTll and with any 
max (IAD > 0.05. 



B. Enhancement factor function of temperature 

Tables IXTVl to IXlXl are the same as Tables IVjtolXll 
except that the enhancement factor / is a function of 
temperature, given by Eq. 21 

A general result that applies to all the rates shown in 
Tables IXIVI to IXIXI is that the sensitivities drop (in ab- 
solute value) strongly as compared to the case of fixed 
enhancement factor. Very few reactions have \Di\ > 0.3 
when we compute the enhancement factor using Eq. |4] 
The list of species sensitive to the rates of (n, 7) reac- 
tions is much shorter, and the only product species with 
\Di\ > 0.3 in this list is 21 Ne. Radiative captures of pro- 
tons suffer as well from a reduction of their influence on 
the supernova yields: 30 Si + p ^ 31 P + 7 continues being 



pernovae. The main reason is that, in our models, 44 Ti 
is made in moderately neutronized matter (77 > 0.001). 
In QSE, while the composition of the Si-group is nearly 
independent of the neutron excess, that of the Fe-group 
is strongly affected 53]. Thus, an increase in the neutron 
excess favours the equilibrium linking of 44 Ti (a r\ = 
nucleus) to the Si-group, leading to a low sensitivity of 
its abundance to the rate of the reactions within the Fe- 
group. Such a progressive decrease of the importance of 
the Fe-group reactions for the production of 44 Ti as 77 in- 
creases, can also be deduced from comparison of Tables 
4, 7, and 8 in 



TABLE XII. Rates that influence the yields of each product species, from Carbon to ChlorineQ 



Nucleus 


End product^ 


Production 
factor c 


Reaction and Di 


ib 




9.49 


2U Ne(«, 7 ):-0.03; 12 C (n, 7 ):+0.02; 12 C (a, 7 ):-0.02 


20 Ne 




1.17 


20 Ne(a,7):-0.29; ie O (a,7):-0.11; 20 Ne(a,p):+0.03; 23 Na(a,p):-0.03; 
24 Mg(a,7):+0.02; 16 (n,7):+0.02; 23 Na(p,7):-0.01 


23 Na 




0.29 


20 Ne(a,p):-0.46; 23 Na(a,p):-0.23; 30 Si(p,7):+0.15; 20 Ne(a,7):+0.12; 27 Al(p,7):+0.08; 
26 Mg(p,7):+0.07; 16 (a,7):+0.07; 24 Mg(Q,7):+0.06; 27 Al(a,p):-0.05; 12 C (n, 7 ):+0.04 


24 Mg 




22.5 


20 Ne(a,7):+0.70; 24 Mg(a,7):-0.42; 24 Mg(a,p):-0.17; 27 Al(p,7):-0.11; 30 Si(p,7):-0.09; 
27 Al(a,p):-0.08; ie O (a,7):-0.07; 24 Mg(n,7):-0.05; 12 C (n, 7 ):-0.05; 26 Mg(a,n):-0.04 


25 Mg 




0.14 


20 Ne(a,7):+0.37; 2j Mg(p,7):-0.29; 2j Mg(a,n):-0.22; 24 Mg(n,7):+0.22; 30 Si(p,7):+0.20; 
24 Mg(a,p):-0.17; 25 Mg(n,7):-0.17; 27 Al(a,p):-0.16; 29 Si(a,n):+0.15; 26 Mg(a,n):+0.10 


26 Mg 




0.22 


26 Mg(p, 7 ):-0.34; 26 Mg(Q,n):-0.34; 23 Na(a,p):+0.31; 30 Si(p,7):+0.23; 20 Ne(a,p):-0.17; 








27 Al(a,p):-0.13; 20 Ne(a,7):+0.12; 27 Al(p,7):+0.10; 24 Mg(a,p):-0.07; 28 Si(a,p):+0.06 
26 Al(p, 7 ):-0.67; 24 Mg(n,7):+0.46; 25 Mg(a,n):-0.39; 25 Mg(p,7):+0.38; 20 Ne(a,7):+0.37; 
30 Si(p,7):+0.26; 29 Si(Q,n):+0.21; 24 Mg(a,p):-0.19; 27 Al(a,p):-0.17; 26 Mg(a,n):+0.15 


26 Al 


26 Mg (7.2 x 10 5 yr) 


9.4 x 10" 4 


27 Al 




4.98 


20 Ne(a,7):+0.48; 27 Al(a,p):-0.44; 27 Al(p,7):-0.35; 24 Mg(a,7):-0.28; 24 Mg(Q,p):+0.26; 
30 Si(p,7):+0.09; 27 Al(a,n):-0.06; 16 (a, 7 ):-0.06; 26 Mg(a,n):-0.05; 28 Si(a,p):+0.04 


28 Si 




220. 


13 N (q, P ):-0.09; 20 Ne(a,7):+0.07; 12 C (a,7):-0.06; 23 Na(a,p):+0.03; 
16 (a, 7 ):-0.02; 12 C (n, 7 ):-0.01; 28 Si(a,7):-0.01 


29 Si 




8.26 


20 Ne(a,7):+0.21; 27 Al(a,n):+0.18; 29 Si(a,n):-0.16; 30 Si(p,7):+0.13; 27 Al(a,p):-0.11; 
23 Na(a,p):+0.11; 30 Si(p,n):-0.09; 26 Mg(a,n):+0.07; 32 S (n,7):-0.06; 30 Si(a,n):-0.06 


30 Si 




22.2 


20 Ne(a,7):+0.37; 30 Si(p,7):-0.27; 30 Si(a,7):-0.16; 27 Al(a,p):+0.14; 28 Si(a,7):-0.11; 
31 P (a,p):-0.09; 24 Mg(a,7):-0.08; 24 Mg(a,p):+0.08; 12 C (a, 7 ):-0.06; 26 Mg(a,n):-0.05 


31p 




30.0 


28 Si(a,p):-0.20; 30 Si(p,7):+0.19; 28 Si(a,7):-0.10; 30 Si(a,7):-0.09; 31 P (a,p):-0.06; 
24 Mg(a,7):+0.06; 27 Al(p, 7 ):+0.04; 27 Al(a,p):+0.04; 31 P (p,7):-0.03; 28 Si(n, 7 ):-0.03 


32 g 




240. 


13 N (q, P ):-0.04; 20 Ne(a,7):-0.03; 28 Si(a,7):+0.01; 16 (a,7):-0.01 


33g 




35.4 


32 S (n, 7 ):+0.18; 30 Si(a,n):-0.17; 20 Ne(a,7):-0.17; 30 Si(p,7):-0.16; 27 Al(a,p):+0.12; 
28 Si(a,7):+0.09; 24 Mg(n,7):+0.07; 23 Na(a,p):-0.06; 28 Si(n,7):+0.05; 27 Al(a,n):-0.05 


34 g 




36.5 


30 Si(p,7):+0.14; 20 Ne(a,7):-0.10; 30 Si(a, 7 ):+0.09; 24 Mg(a,p):-0.09; 28 Si(Q,7):+0.07; 
31 P (a,p):+0.07; 27 Al(p, 7 ):+0.06; 29 Si(a,n):+0.05; 24 Mg(a,7):+0.05; 32 S (a,p):+0.05 


35 CI 




11.0 


20 Ne(a,7):-0.35; 30 Si(a,7):+0.20; 31 P (a,p):+0.12; 30 Si(p,7):+0.11; 24 Mg(a,p):-0.07; 
24 Mg(a,7):+0.07; 32 S (a,p):-0.07; 29 Si(a,n):+0.07; 27 Al(p,7):+0.07; 28 Si(Q,7):+0.05 



a Only the species with mass fraction greater than 10~ 5 or that are interesting radioactive isotopes and with max(|Di|) > 0.01 are shown here. For each nucleus, we show a 
maximum of ten reactions. 

b For radioactive nuclei, here it is shown the end product of the disintegration chain as well as the longest half-life in the decay chain. 

c Mass fraction of the species in the supernova ejecta normalized to its solar mass fraction. For radioactive nuclides the normalization is with respect to the solar mass fraction of 
the end product of the disintegration chain. 
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FIG. 12. (Color online) Chart of the most influential reactions compared to the maximum mass flow path. Left: Reactions 
with max (|-Dj|) > 0.3 in Tables IVTl to IXII The type of each reaction can be deduced from the differences of atomic and baryonic 
number of the nuclei it connects or from the color in which they are drawn: (n, 7) reactions in red, (p, 7) reactions in cyan, 
(p, n) reactions in magenta, (a, 7) reactions in blue, (a, n) reactions in black, and (a, p) reactions in green. We recall that 
direct and inverse reaction rates were modified simultaneously. Right: Reactions that carry a large mass flow. In this plot, the 
mass flow has been color coded according to the value of dJtjk' red for DJljk > 0.5 Mq, magenta for 0.5 > ZDljk > 0.4 Mq, blue 
for 0.4 > Ttjk > 0.3 Mq, green for 0.3 > m jk > 0.2 Mq, and cyan for 0.2 > m jk > 0.01 Mq. The chart of most influential 
reactions has been superimposed as black thin lines. 



TABLE XIII. Rates that influence the yields of each product species, from Argon to Nickel (continuation of Table IXTlj) Fl 



Nucleus 


End product^ 


Production 
factor c 


Reaction and Di 


3b Ar 




270. 


2U Ne(a, 7 ):-0.09; i2 C (a, 7 ):+0.04; 13 N (a,p):+0.03; 
23 Na(a,p):-0.03; 44 Ti(a,p):-0.02; 23 Na(p, 7 ):-0.02 


37 Ar 


37 CI (35 d) 


7.37 


20 Ne(a, 7 ):-0.25; 34 S (a,n):-0.15; 36 Ar(n, 7 ):+0.11; 12 C (a, 7 ):+0.07; 30 Si(a, 7 ):+0.06; 
35 Cl(p, 7 ):+0.06; 33 S (a, 7 ):+0.06; 13 N (a,p):-0.05; 24 Mg(a, 7 ):+0.05; 42 Ca(a, 7 ):-0.05 


38 Ar 




13.6 


20 Ne(a, 7 ):-0.22; 35 Cl(a,p):+0.11; 34 S (q, 7 ):+0.09; 30 Si(p, 7 ):+0.09; 30 Si(a, 7 ):+0.08; 
24 Mg(a, 7 ):+0.08; 28 Si(a, 7 ):+0.06; 13 N (a,p):-0.05; 42 Ca(a, 7 ):-0.05; 24 Mg(a,p):-0.05 


39 K 




7.81 


20 Ne(a, 7 ):-0.35; 12 C (a, 7 ):+0.16; 35 Cl(a,p):+0.09; 42 Ca(a, 7 ):-0.07; 24 Mg(a, 7 ):+0.07; 
30 Si(a, 7 ):+0.07; 39 K (p, 7 ):-0.07; 34 S (q, 7 ):+0.05; 45 Sc(p, 7 ):-0.05; 30 Si(p, 7 ):+0.04 


40 Ca 




370. 


20 Ne(a, 7 ):-0.15; 13 N (a,p):+0.11; 12 C (a, 7 ):+0.10; 23 Na(a,p):-0.05; 23 Na(p, 7 ):-0.03; 
16 (a, 7 ):+0.03; 44 Ti(a,p):-0.03; 21 Na(a,p):-0.01; 12 C (n, 7 ):+0.01 


44 Ti 


44 Ca (60 yr) 


14.3 


20 Ne(a, 7 ):-0.17; 44 Ti(a,p):-0.17; 13 N (a,p):+0.16; 12 C (a, 7 ):+0.14; 40 Ca(a, 7 ):-0.08; 
23 Na(a,p):-0.06; ie O (a, 7 ):+0.04; 23 Na(p, 7 ):-0.03; 21 Na(a,p):-0.02; 12 C (n, 7 ):+0.01 


48 y 


48 Ti (16 d) 


96.8 


13 N (a,p):+0.11; 4S Cr(a,p):-0.10; 20 Ne(a, 7 ):-0.09; 44 Ti(a,p):+0.06; 12 C (a, 7 ):+0.05; 
23 Na(a,p):-0.04; ie O (a, 7 ):+0.03; 49 V (p,n):-0.02; 52 Fe(a,p):-0.01; 23 Na(p, 7 ):-0.01 


49y 


49 Ti (330 d) 


59.8 


20 Ne(a, 7 ):-0.08; 49 V (p,n):-0.08; 48 Cr(a,p):-0.07; 13 N (a,p):+0.07; 44 Ti(a,p):+0.06; 
12 C (a, 7 ):+0.05; 53 Mn(p,n):-0.04; 23 Na(a,p):-0.02; 49 V (p, 7 ):-0.02; ie O (a, 7 ):+0.02 


50 Cr 




100. 


51 Mn(p, 7 ):-0.11; 55 Co(p, 7 ):+0.10; 13 N (a,p):-0.06; 20 Ne(a, 7 ):+0.04; 53 Mn(p,n):-0.03; 
23 Na(a,p):+0.02; 52 Fe(a,p):-0.02; 50 Cr(n, 7 ):-0.02; 52 Mn(p, 7 ):-0.02; 20 Ne(a,p):-0.01 








51 Cr 


51 V (28 d) 


73.8 


51 Mn(p, 7 ):-0.17; 20 Ne(a, 7 ):-0.05; 13 N (a,p):+0.03; 50 Cr(p, 7 ):+0.03; 53 Mn(p,n):-0.03; 
44 Ti(a,p):+0.03; 52 Fe(a,p):-0.02; 12 C (a, 7 ):+0.02; ie O (a, 7 ):+0.01 


52 Mn 


52 Cr (5.6 d) 


280. 


13 N (a,p):+0.12; 20 Ne(a, 7 ):-0.07; 12 C (a, 7 ):+0.05; 16 (a, 7 ):+0.04; 23 Na(a,p):-0.04; 
44 Ti(a,p):+0.03; 52 Fe(a,p):-0.03; 53 Mn(p,n):-0.03; 42 Ca(a, 7 ):+0.01; 45 Sc(p, 7 ):+0.01 


53 Mn 


53 Cr (3.7 x 10 6 yr) 


240. 


53 Mn(p,n):-0.12; 13 N (a,p):+0.07; 20 Ne(a, 7 ):-0.05; 12 C (a, 7 ):+0.03; 52 Mn(p,n):+0.03; 
ie O (a, 7 ):+0.02; 23 Na(a,p):-0.02; 44 Ti(a,p):+0.02; 52 Fe(a,p):-0.02; 53 Mn(p, 7 ):-0.01 


55 Fe 


55 Mn (2 7 vrl 


260. 


55 Co(p, 7 ):-0.11; 13 N (a,p):+0.04; 20 Ne(a, 7 ):-0.03; 
12 C (a, 7 ):+0.02; ie O (a, 7 ):+0.02; 23 Na(a,p):-0.01 


57 Co 


57 Fe (270 d) 


300. 


57 Co(p,n):-0.07; 13 N (q, P ):+0.01 


56 Ni 


56 Fe (77 d) 


360. 


13 N (a,p):+0.03; 20 Ne(a, 7 ):-0.01; 12 C(a, 7 ):+0.01 


58 Ni 




410. 


57 Co(p,n):+0.01 


59 Ni 


59 Co (7.6 x 10 4 yr) 


83.4 


59 Cu(p, 7 ):-0.18; 56 Ni(a,p):+0.13; 58 Ni(n, 7 ):-0.01; 57 Co(p,n):+0.01 


60 Ni 




210. 


56 Ni(a,p):+0.13; 59 Cu(p, 7 ):+0.07; 56 Ni(a, 7 ):+0.01 


61 Ni 




110. 


57 Ni(a,p):+0.16; 57 Co(p,n):-0.08; 57 Ni(a, 7 ):+0.07; 60 Cu(p, 7 ):+0.02; 21 Na(a,p):-0.01 


62 Ni 




160. 


58 Ni(a, 7 ):+0.12; 58 Ni(a,p):+0.10; 60 Zn(a,p):+0.05; 57 Co(p,n):+0.04; 
21 Na(a,p):-0.02; 59 Cu(a,p):+0.02; 61 Cu(p, 7 ):+0.01 



a Only the species with mass fraction greater than 10~ 5 or that are interesting radioactive isotopes and with max(|Di|) > 0.01 are shown here. For each nucleus, we show a 
maximum of ten reactions. 

For radioactive nuclei here it is shown the end product of the disintegration chain as well as the longest half-life in the decay chain. 
c Mass fraction of the species in the supernova ejecta normalized to its solar mass fraction. For radioactive nuclides the normalization is with respect to the solar mass fraction of 
the end product of the disintegration chain. 
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TABLE XIV. Sensitivity of the nucleosynthesis to the rate of (n, 7) reactions with the parent nuclide given in the first column, 
with enhancement factor given by Eq. flFl 



"PaTPTTt" tiiipIiHp 

± Cll Cllu 11 LLL^llvlty 


Nuclei with ID 1 > S 

J. 1 Li l^lCl Vv lull \JSi ^ U 1 iJ 


Nuclei with 3 > ID- > 05 

J. N LI ^ 1 \1' 1 Vv lull U 1 iJ J_^'2 UiUy 


01 






55 Fe 




55 Mn 
1V1I1 


32 g 




32,33p 33,35g 


36 Ar 




37 Cl, 37 Ar 


24 Mg 




21 Ne, 25 Mg, 26 Al, 35 S 


25 Mg 




17 0, 21 Ne, 25 Mg, 26 Al 


56 Fe 




57 Fe 


16 




17 


46 Ti 




46 Ti 


20 Ne 


21 Ne 




!2 C 




21 Ne 


31p 




32 p 



a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table ITTljl and with any 
max (| |) > 0.05. 



TABLE XV. Sensitivity of the nucleosynthesis to the rate of (p, 7) reactions with the parent nuclide given in the first column, 
with enhancement factor given by Eq. HFI 



Parent nuclide 


Nuclei with \D t \ > 0.3 


Nuclei with 0.3 > \D t \ > 0.05 


29 Si 




2b Al, 3b S, 43 Ca 


35 CI 




37 C1 


30 Si 


35 S 


23 Na, 25 ' 26 Mg, 26 Al, 29 ' 30 Si, 31 ' 32 P, 33 ' 34 S, 35 Cl, 43 Ca, 47 Ti 


27 Al 




23 Na, 24 - 26 Mg, 26 ' 27 Al, 32 P, 35 S, 43 Ca, 47 Ti 


59 Cu 




59 Ni 


56 Co 




56 Co 


33g 




43 Ca 


51 Mn 




50,51 Cr 


47y 




47 Ti 


26 Al 


26 Al 




59 Co 




59 Co 


26 Mg 




26 Mg, 26 Al, 35 S, 43 Ca, 47 Ti 


25 Mg 




21 Ne, 25 Mg, 26 Al, 35 S 


58 Fe 




58 Fe 


45 Sc 




45 Sc 


B2 Cu 




63 Cu, 65 Zn 


37 C1 




37 C1 


17 F 




14 N 


64 Ga 




63 Cu, 65 Zn 



a The reactions listed are those that processed more than 10 6 in the reference model (see Table ITTljl and with any 
maxflDil) > 0.05. 



TABLE XVI. Sensitivity of the nucleosynthesis to the rate of (p, n) reactions with the parent nuclide given in the first column, 
with enhancement factor given by Eq. l4f*l 



Parent nuclide 


Nuclei with |A| > 0.05 


Parent nuclide 


Nuclei with |A > 0.05 


M Mn 


b3 Mn 


3^p 


3^p 


30 Si 


35 g 


45 Sc 


45 Sc 


58 Co 


58 Co 


37 C1 


37 C1 


54 Mn 


54 Mn 


62 Cu 


63 Cu, 65 Zn 


55 Mn 


55 Mn 


47 Ti 


46 Ti 



a The reactions listed are those that processed more than 10 6 Mq in the reference model (see Table ITTTll and with any 
maxflDil) > 0.05. 
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TABLE XVII. Sensitivity of the nucleosynthesis to the rate of (a, 7) reactions with the parent nuclide given in the first column, 
with enhancement factor given by Eq. [W\ 



Parent nuclide 


Nuclei with |A| > 0.3 


Nuclei with 0.3 > |A| > 0.05 


28 Si 




3U Si 


32 g 




37 C1 


20 Ne 




20 < 21 Ne, 23 Na, 24 - 26 Mg, 26 ' 27 Al, 29 ' 30 Si, 32 ^^ 45 Sc, 

44-47rpj 






21 Ne, 23 Na 


Mg 


43 s < 
Ca 


24n iT 27 a i 35n 47m. 

Mg, Al, S, Ti 


58 Ni 




62 Ni, 63 Cu, 64 Zn 


12 C 




39 K, 41 ' 42 Ca, 45 Sc, 44 ' 46 Ti 


33 g 




37 C1 


30 Si 




30gj 32,33p 34,35g 35q 


41 Ca 




43 Ca, 47 Ti ' 


42 Ca 




46 Ti 


62 Zn 




66 Zn 



a The reactions listed are those that processed more than 10 6 M0 in the reference model (see Table iHlll and with any 
max(|Di|) > 0.05. 



TABLE XVIII. Sensitivity of the nucleosynthesis to the rate of (a, n) reactions with the parent nuclide given in the first column, 
with enhancement factor given by Eq. [4n 



Parent nuclide 


Nuclei with |A| > 0.3 


Nuclei with 0.3 > \D t \ > 0.05 


29 Si 




2b Mg, 2b Al, 29 Si, 3b S, 43 Ca 


33 g 




37 Cl, 43 Ca 


27 A1 


43 Ca 


25 Mg, 29 Si, 32 ' 33 P, 35 S, 47 Ti 


30 Si 




33 P, 33 - 35 S, 37 Cl, 43 Ca, 47 Ti 


25 Mg 




17 0, 21 Ne, 25 Mg, 26 Al, 35 S 


34g 




37 C1 37 Ar 


26 Mg 




21 Ne', 25 ' 26 Mg, 26 Al, 32 P, 35 S, 43 Ca, 47 Ti 


41 Ca 




43 Ca, 47 Ti 


21 Ne 


21 Ne 




17 




17 0, 21 Ne 



a The reactions listed are those that processed more than 10 6 M0 in the reference model (see Table [TTTT) and with any 
maxflDil) > 0.05. 



TABLE XIX. Sensitivity of the nucleosynthesis to the rate of (a, p) reactions with the parent nuclide given in the first column, 
with enhancement factor given by Eq. l4pl 



Parent nuclide 


Nuclei with \Di\ > 0.05 


Parent nuclide 


Nuclei with |A| > 0.05 


28 Si 


31,33p 


44 Ti 




30p 


43 Ca 


57 Ni 


61 Ni, 65 Zn 


32 g 


32,33p 35g ^^Q]^ 


23 Na 


21 Ne, 23 Na, 26 Mg, 29 Si 43 Ca, 47 Ti 


27 Al 


2S,26^ g) 26,27 A1) 29,30 Si) 32 Pj 33,35 S) 43 Ca 


41 Ca 


43 Ca 


31p 


33 P, 35 S, 35 C1 


34g 


37 CI 


24 Mg 


24,25 Mg] 26,27 A1] 35 S) 43 Ca] 47 Ti 


39 Ca 


43 Ca 


56 Ni 


43 Ca, 59 ' 60 Ni, 63 Cu, 64 Zn 


30 Si 


33p 


39 K 


43 Ca 


42 Ca 


46 Ti 


40 Ca 


43 Ca 


62 Zn 


66 Zn 


29 Si 


32 p 


61 Zn 


65 Zn 


13 N 


45 Sc, 44 Ti 


17 Ne 


14 N 


35 CI 


38 Ar 


43 Ti 


43 Ca, 47 Ti 


20 Ne 


17 0, 21 Ne, 23 Na, 26 Mg 


47 Cr 


47 Ti 


58 Ni 


62 Ni, 63 Cu, 64 Zn 







a The reactions listed are those that processed more than 10 6 M0 in the reference model (see Table ITTT)| and with any 
maxflDil) > 0.05. 
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the most influential reaction, but the number and im- 
portance of the species affected by its rate is much lower 
than with a fixed rate enhancement factor. The only 
species with |A| > 0.3 are 26 Al and 35 S. The influence of 
(p, n) reactions on the supernova yields is marginal when 
Eq. [4] is used to determine the enhancement factor of the 
rates, and the same applies to (a, p) reactions. Among 
the (a, 7) reactions, the capture on 20 Ne continues being 
the reaction with the largest list of product species with 
\Di\ > 0.05. The only species with |_D,-| > 0.3 due to vari- 
ations on the rate of this type of reactions is the trace 
species 43 Ca. Finally, the only species with \Di\ > 0.3 
with respect to variations of (a, n) reaction rates in Ta- 
ble KVIIII are 21 Ne and 43 Ca. 

The rate enhancement factor computed with Eq.Q] dif- 
fers most from the fixed enhancement factor at temper- 
atures T > 3-5 x 10 9 K. Thus, it has a stronger effect 
on reactions whose main role is played at high temper- 
atures, as highlighted by comparing, once more, the list 
of reactions (parent nuclei) in Tables IXIVI to IXIXI with 
that in Tables IVfl to 1X11 For instance, as most (n, 7) 
reactions influence the yields at temperatures of order 
2 x 10 9 < T < 4 x 10 9 K, the list of reactions in Ta- 
bles [VI] and [XlV] is quite similar. On the other side, the 
list of (p,7) reactions in Table IXVl is much shorter than 
in Table IVffl because the reduction of the enhancement 
factor at high temperatures affects most reactions with 
Fe-group nuclei, while reactions with IMEs are less af- 
fected. 

Tables |XX] and |XXT] are similar to Tables [XII] and KnU 

except that the enhancement factor / is now a function 
of temperature, given by Eq. [4] The maximum \Di\ 
achieved with Eq. U for a given species is in general a 
factor of two smaller than when using a fixed rate en- 
hancement factor, while all the species with a production 
factor larger than 100 have maximum \Di\ < 0.1. 



C. Sensitivity to different rate prescriptions 

In this Section, we analyze the changes in the yields 
obtained using different prescriptions for the rates of a 
few selected reactions. To this end, we have accessed the 
JINA REACLIB Database to compare the most recent 
rates for each one of the selected reactions. We discuss 
in the following the prescriptions for the reactions that 
appear in Tables IVT1 to IXll with maximum \Di\ > 0.3. 
The results are presented in Table IXXIII in the form of 
percent variations of the yield of product species when 
two different prescriptions are used for each reaction rate. 
We give in the table as well the sources of the rates of 
each reaction. The reference rate (i.e. that used in the 
denominator of the calculation of the relative variation 
of the yield) is always that cited in second place in the 
table. All the references that appear in this Section are 
taken from the JINA webpage. 



1. (n, 7) reactions 

The rates of the reactions we consider are fits to either 
theoretical or experimental results published in [42], [56l — 

The discrepancy between the different rates of the 
reactions 32 S + n ^ 33 S + 7, 28 Si + n ^ 29 Si + 7, 
24 Mg + n ^ 25 Mg + 7, 25 Mg + n ^ 26 Mg + 7, and 
20 Ne + n ^ 21 Ne + 7 computed from the above refer- 
ences is less than a factor of ten for T > 10 9 K. As this 
uncertainty is within the range explored in Section [IV Al 
we do not deem it necessary to further analyze these re- 
action rates. 

On the other hand, the rate of the reaction 
16 + n ?± 17 + 7 computed from the two references in 
JINA, [56| and [58| . shows a discrepancy of more than 
two orders of magnitude between these two cases. We 
have computed the nucleosynthesis of our SNIa model 
with both rates and compared the results in the first 
row of Table IXXIII Aside from the trace product 17 0, 
whose yield decreases by two orders of magnitude when 
using the rate from [58| . the effect on each abundance 
is smaller than 27%. We conclude that (n,j) reaction 
rates, in general, are not critical for obtaining accurate 
yields from SNIa models. 

2. (p, 7) reactions 

The rates of the reactions we consider are fits to either 
theoretical or experimental results published in [42], l57l . 

The discrepancy between the different rates of the re- 
actions 25 Mg + p <± 26 Al + 7, 26 Mg + p ^ 27 Al + 7. 
26 Al + p <=i 27 Si + 7, 27 A1 + p <=> 28 Si + 7, and 30 Si + p «=* 
31 P + 7 computed from the above references is less than 
a factor of ten for T > 10 9 K, well within the range ex- 
plored in Section ITVAl The reaction 45 Sc + p <=* 46 Ti + 7, 
which contributes to the linking of QSE groups in silicon 
burning, is only evaluated in JINA through three some- 
what different theoretical models. In the temperature 
range of interest, the rates given by these models match 
perfectly. 

The rate of the reaction 17 F + p <=2 18 Ne + 7 is given in 
JINA for two different fits to experimental rates, which 
differ by more than one order of magnitude for temper- 
atures in the range 10 9 < T < 10 10 K. We have recom- 
puted the nucleosynthesis with both evaluations of the 
rate of this reaction and show the results in the second 
row of Table IXXIII Only the yield of 14 N, a marginal 
product of the supernova nucleosynthesis, changes by 
more than 10%. 

The reaction 30 Si+p <=^ 31 P+7 has the largest D.;, con- 
sequently we have recomputed the nucleosynthesis using 
the three most recent evaluations of its rate from JINA. 
The results are shown in the third and fourth rows of 
Table IXXIII in which we have taken as a reference the 
recommended rate from (60| . which is compared to two 



TABLE XX. Rates that influence the yields of each product species, from Carbon to Chlorine, when the enhancement factor is given by Eq. |4"F1 



Nucleus 


Reaction and Di 


IbQ 


2U Ne(a,7):-0.01; 12 C (n,7):+0.01 


20 Ne 


20 Ne(a,7):-0.15; 23 Na(a,p):-0.01; 16 (a,7):-0.Ql; 24 Mg(Q, 7 ):+0.01; ie O (n,7):+0.01 


23 Na 


20 Ne(a,p):-0.24; 23 Na(a,p):-0.13; 20 Ne(a,7):+0.08; 30 Si(p,7):+0.07; ls O (a,7):+0.06; 27 Al(p,7):+0.05; 24 Mg(a, 7 ):+0.04; 26 Mg(p,7):+0.04; 




27 Al(a,p):-0.03; 12 C (n,7):+0.02 


24 Mg 


20 Ne(a,7):+0.25; 24 Mg(a, 7 ):-0.18; 24 Mg(a,p):-0.09; 27 Al(p,7):-0.06; 30 Si(p, 7 ):-0.05; 27 Al(a,p):-0.04; 16 (a, 7 ):-0.03; 12 C (n, 7 ):-0.03; 




24 Mg(n,7):-0.02; 29 Si(a,n):-0.02 


25 Mg 


20 Ne(a,7):+0.19; 25 Mg(p, 7 ):-0.17; 25 Mg(a,n):-0.13; 30 Si(p,7):+0.12; 24 Mg(n,7):+0.10; 29 Si(a,n):+0.09; 27 Al(a,p):-0.09; 25 Mg(n,7):-0.09; 




24 Mg(a,p):-0.07; 27 Al(p, 7 ):+0.05 


26 Mg 


26 Mg(p, 7 ):-0.20; 26 Mg(a,n):-0.20; 23 Na(a,p):+0.17; 30 Si(p,7):+0.10; 20 Ne(a,7):+0.10; 27 Al(a,p):-0.07; 27 Al(p, 7 ):+0.07; 20 Ne(a,p):-0.06; 




24 Mg(a,p):-0.03; 24 Mg(a, 7 ):+0.03 


26 Al 


26 Al(p, 7 ):-0.33; 25 Mg(a,n):-0.22; 24 Mg(n, 7 ):+0.18; 2o Ne(a, 7 ):+0.18; 25 Mg(p, 7 ):+0.17; 30 Si(p, 7 ):+0.14; 29 Si(a,n):+0.12; 27 Al(a,p):-0.10; 




24 Mg(a,p):-0.08; 26 Mg(p, 7 ):-0.08 


27 Al 


27 Al(a,p):-0.23; 27 Al(p,7):-0.18; 20 Ne(a,7):+0.17; 24 Mg(a,7):-0.13; 24 Mg(a,p):+0.06; 26 Mg(a,n):-0.03; 27 Al(a,n):-0.03; 16 (a,7):-0.02; 




26 Mg(p,7):+0.02; 30 Si(a,7):-0.02 


28 Si 


13 N (a,p):-0.03; 20 Ne(a,7):+0.02; 12 C (a, 7 ):-0.02; 23 Na(a,p):+0.01 


29 Si 


27 Al(a,n):+0.09; 29 Si(a,n):-0.08; 20 Ne(a,7):+0.07; 30 Si(p, 7 ):+0.06; 27 Al(a,p):-0.06; 23 Na(a,p):+0.06; 30 Si(p,n):-0.04; 26 Mg(a,n):+0.04; 




26 Mg(p, 7 ) :-0.03; 32 S (n, 7 ):-0.03 


30 Si 


30 Si(p,7):-0.14; 20 Ne(a,7):+0.13; 30 Si(a,7):-0.08; 27 Al(a,p):+0.06; 28 Si(a, 7 ):-0.05; 31 P (a,p):-0.04; 32 S (a,p):-0.03; 12 C (a, 7 ):-0.02; 




26 Mg(a,n):-0.02; 24 Mg(a, 7 ):-0.02 


31p 


30 Si(p,7):+0.07; 28 Si(a,p):-0.06; 28 Si(a,7):-0.05; 20 Ne(a, 7 ):-0.04; 30 Si(a, 7 ):-0.04; 31 P (a,p):-0.03; 24 Mg(a, 7 ):+0.03; 27 Al(p,7):+0.02; 




27 Al(a,p):+0.02; 32 S (a,p):-0.02 


32 s 


13 N (a,p):-0.02; 20 Ne(a,7):-0.01 


33 s 


20 Ne(a,7):-0.11; 32 S (n,7):+0.09; 30 Si(a,n):-0.08; 27 Al(a,p):+0.06; 30 Si(p, 7 ):-0.06; 28 Si(a,7):+0.05; 24 Mg(n,7):+0.04; 23 Na(a,p):-0.03; 




27 Al(a,n):-0.03; 28 Si(n, 7 ):+0.03 


34 s 


20 Ne(a,7):-0.12; 30 Si(p, 7 ):+0.05; 30 Si(a,7):+0.05; 28 Si(a,7):+0.04; 31 P (a,p):+0.04; 32 S (a,p):+0.03; 27 Al(p,7):+0.03; 24 Mg(a,p):-0.03; 




29 Si(a,n):+0.03; 24 Mg(a,7):+0.02 


35 CI 


20 Ne(a,7):-0.15; 30 Si(a,7):+0.11; 31 P (a,p):+0.06; 30 Si(p,7):+0.05; 27 Al(p,7):+0.04; 29 Si(a,n):+0.04; 24 Mg(a, 7 ):+0.02; 28 Si(a,7):+0.02; 23 Na(a,p):-0.02; 




24 Mg(a,p):-0.02 



a Only the species with mass fraction greater than 10 5 or that are interesting radioactive isotopes and with max(|Di|) > 0.01 are shown here. For each nucleus, we show a 
maximum of ten reactions. 
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TABLE XXI. Rates that influence the yields of each product species, from Argon to Nickel, when the enhancement factor is given by Eq. [4] (continuation of Table |XX)) F1 



Reaction and D, 
2U Ne(a,7):-0.04; 
20 Ne(a,7):-0.09; 34 S (a,n, 
27 Al(a,p):+0.02; 37 Cl(p,n 

20 



i2 C (a,7):+0.02; 13 N (q, P ):+0.01 
34 S (a,n):-0.08; 36 Ar(n, 7 ):+0.06; 12 C (a, 7 ) 
0.02 
-0.05 



-0.04: 



Ne(a,7):-0.10; 35 Cl(a,p) 
13 N (q, P ):-0.02; 29 Si(Q,n):+0.02 
20 Ne(a,7):-0.11; 12 C (a,7):+0.08 
34 S (a, 7 ):+0.02; 30 Si(p, 7 ):+0.02 
20 Ne(a,7):-0.07; 13 N (a,p):+0.04 
20 Ne(a,7):-0.07; 44 Ti(a,p) 
48 Cr(a,p):-0.04; 13 N (a,p):+0.04 
49 V (p,n):-0.04; 20 Ne(a, 7 ):-0.03; 
51 Mn(p, 7 ):-0.05; 55 Co(p,7):+0.03 



34 S (a, 7 ):+0.05; 

; 35 Cl(a,p):+0.04; 

; 12 C (a, 7 ):+0.04 
0.07; 13 N (a,p):+0.07; 
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30 Si(p,7):+0.04; 

30 Si(a,7):+0.03; 

; 23 Na(a,p):-0.02; 
12 C (a, 7 ):+0.05; 

44 r 



13 N (q, P ):-0.03; 35 Cl(p,7):+0.03; 30 Si(Q,7):+0.03; 33 S (a, 7 ):+0.03; 
30 Si(a,7):+0.04; 24 Mg(a, 7 ):+0.03; 28 Si(a,7):+0.03; 27 Al(p, 7 ):+0.02; 



13 t- 



N (a,p):-0.03; 39 K (p, 7 ):-0.03; 24 Mg(a,7):+0.03; 42 Ca(Q, 7 ):-0.03 



44 Ti(a,p):-0.01 

23 Na(a,p):-0.02; 40 Ca(a,7):-0.02 
I2r< I n , „A. i n no. 23 



Ne(a,7):-0.04; 44 Ti(a,p):+0.03; 12 C (a, 7 ):+0.02 



51 Mn(p, 7 ) :-0.09; 20 Ne(a,7):-0.02; 

20 Ne(a,7):-0.03; 
20 Ne(a, 7 ):-0.02; 



48 Cr(a,p):-0.03; 13 N (a,p):+0.02; 44 Ti(a,p):+0.02; 53 Mn(p,n 
20 Ne(a, 7 ):+0.02; 13 N (a,p):-0.02; 23 Na(a,p):+0.01 



Na(a,p) :-0.02 

0.02; 12 C (q, 7 ):+0.02 



13 N (q, P ):+0.04 
53 Mn(p,n):-0.05; 
55 Co(p, 7 ):-0.05; 
57 Co(p,n) :-0.03 
13 N (q, P ):+0.01 
59 Cu(p, 7 ):-0.09; 
56 Ni(Q,p):+0.07: 
57 Ni(Q,p):+0.09 
58 Ni(Q, 7 ):+0.07 



13 



N (a,p):+0.01: 



55 Co(p,7):-0.01; 
12 C (a,7):+0.02 
13 N (q, P ):+0.02: 
20 Ne(a, 7 ):-0.01 



50 Cr(p,7):+0.01; 52 Fe(a,p):-0.01; 
; 23 Na(a,p):-0.02; 52 Fe(a,p):-0.01: 

12, 



44 Ti(Q,p):+0.01 
44 Ti(Q,p):+0.01 



2 C (a,7):+0.01 



56 Ni(a,p):+0.08 
59 Cu(p, 7 ):+0.02 

57 Ni(a,7):+0.04; 57 Co(p,n):-0.03 
58 Ni(a,p):+0.06; 60 Zn(a,p):+0.02; 



7 Co(p,n):+0.01 



Only the species with mass fraction greater than 10 5 or that are interesting radioactive isotopes and with max(|Di|) > 0.01 are shown here. For each nucleus, we show a 
maximum of ten reactions. 
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TABLE XXII. Sensitivity of the nucleosynthesis to different rate prescriptions for the reactions in the first column. 



Reaction 



5 0+n^ 1Y + 7 ^ 
7 F + p <^ 18 Ne + ~^ 
:, Si + p <=> 31 P + 7- f 

:, Si + p^ 31 P + 7l 



Product nuclei and percent variation^ 

ly O:-99%; 14 N:-27%; 2b Mg:-25%; 22 Ne:+22%; 43 Ca:+18%; 21 Ne: + 15%; 

1 A . ^ 



J Ne + a : 



Ne + a : 



'Mg + 71 



Mg + 7- 



4 Mg + a ^ 27 Al + p^ 



S:+14%: 26 Al:-ll% 



14 N:+79% 

3 
31 



5 S: + I00%; 32 P:+36%; 43 Ca:+33%; 30 Si:-30%; 26 Al:+29%; 25 Mg:+28%; 29 Si:+23%; 26 Mg:+22% 
1 P:+18%; 33 S:-I7%; 34 S: + I6%; 35 C1:+I3%; 23 Na:+I2%; 38 Ar:+10% 

35 S:+68%; 32 P:+26%; 43 Ca:+24%; 30 Si:-22%; 26 Al:+2f%; 25 Mg:+20%; 29 Si:+16%; 26 Mg:+16%; 

s Ti:-74%; 44 Ca:-73%; 26 Al:+72%; 30 Si:+7f%; 



33 S:-I4%; 31 P:+13%; 34 S: + I2% 
24 Mg:+I56%; 27 Al:+90%; 40 K:-80%; 43 Ca:-75% 



35 S:-57%; 36 S:+53%; 3W 25 
45 Sc:-39%; 20 Ne:-38%; 39 K:-38%- 37 
33 S:-30%; 29 Si:+29%; 37 

26at„. i iioy. 28q- i -|n<K. 31 



Cl:-53%; 25 Mg:+52%; 33 P:+49%; 



2 Ca:-49%; 38 Ar:-43%; 47 Ti:-41%; 
Ar:-37%; 41 Ca:-37%; 46 Ti:-36%; 32 P:+34%; 34 S:-34%; 
S:-3()':-r: 29 Si:+29%; 37 Cl:-29%; 44 Ti:-27%; 23 Na:-21%; 14 N:-16%; 21 Ne:-f6%; 40 Ca:-16%; 



Mg:+ll%; 28 Si:+10%; 31 P 

40t^. coOZ. 43 



-10% 



Mg:+84%; 40 K:-58%; 43 Ca:-56% 



7 Al:+52%; 48 Ti:-51%: 



30 Si:+48%; 44 Ca:-48%; 35 S:-37%; 



Ca:-37%; 45 Sc:-37%; 33 P:+36%; 36 S:+35%; 35 Cl:-35%; 47 Ti:-32%; 44 Ti:-30%; 26 Al:+29%, 
P:+29%; 38 Ar:-29%; 41 Ca:-28%; 39 K:-26%; 46 Ti:-26%; 37 Ar:-24%; 25 Mg:+2I%; 29 Si:+21%; 
S:-20%; 40 Ca:-I9%; 20 Ne:-16%; 33 S:-f6%; 14 N:-I5%; 37 C1:-I4%; 28 Si: + I2%; 48 V:-II%; 
S:+52%; 24 Mg:+44%; 26 Al:+35%; 25 Mg:+26%; 27 A1:-19%; 37 C1:+I6%; 43 Ca:+I5%; 34 S:+13%; 
P: + I2%; 26 Mg:+12%; 35 C1:+12%; 32 P:+10% 



a We only show species for which the percent variation is largest than 10% in absolute value and whose yield is > 10 8 Mq. 
b Rate from |58f vs. rate from [56ll 

c The reference rate from [kJ is compared to a rate that incorporates several contributions from M. Wiescher, as given in the JINA 

Database 
d Rate from |4'J vs. rate from 6C | 
c Rate from 57] vs. rate from 6C | 
f Rate from |42| vs. rate from 68 1 
g Rate from | 61l | vs. rate from 68 1 
h Rate from }62l | vs. rate from 6C | 



other evaluations due to [42j and [57J. In the temper- 
ature range of interest, the rates computed from these 
sources differ by less than a factor of three, and the same 
applies to the rate computed from the older reference 
j63| [. The species whose yields are most sensitive to the 
different prescriptions for this rate are more or less the 
same as already noted in Tables I VIII and IXV1 However, 
the changes in the yields are more consistent with those 
shown in Table IXVl indicating that the use of the en- 
hancement factor function of temperature, as in Eq. 21 
might describe better the actual uncertainties in the nu- 
cleosynthesis than using a fixed enhancement factor, at 
least at the current level of knowledge of this reaction 
rate (but see Sections IIV Dl and lv | . 



3. (p, n) reactions 



Only two (p, n) reactions in Table IVIIII have 
any |A| > 0.3, these are 32 P + p^ 32 S + n and 
37 C1 + p T± 37 Ar + n. In the JINA library there are only 
theoretical rates of these reactions, all of them obtained 
from the NON-SMOKER code, using different nuclear in- 
puts. In the range of temperatures of interest, the differ- 
ent rates for these reactions match each other perfectly. 



:!() 



4- (a, 7) reactions 
The rates of the reactions 24 Mg + a 



28c 



7 and 

Si + a <=2 34 S + 7 are derived from different evaluations 
obtained with the NON-SMOKER code with different 
nuclear inputs, and from [63j. All these rates, for a given 
reaction, agree within a factor smaller than the enhance- 
ment factor we have explored earlier in this paper, hence 
we do not continue with the analysis of these two reac- 
tions. 

The different prescriptions for the rate of the reaction 
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Ne- 
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Mg+7 are from [42J, |5JJ, |61|, |63|, |68| . Within 
the temperature range of interest, these rates show dis- 
crepancies of nearly an order of magnitude. Thus, as this 
reaction is one of the most influential for SNIa nucleosyn- 
thesis, we have recomputed the yields for the three most 
recent prescriptions of its rate. The results are shown 
in the fifth and sixth rows in Table IXXIII The list of 
species whose yields are most sensitive to the prescrip- 
tion used for this reaction is quite similar to that found in 
Table HXl The first point to note is the long list of species 
whose yield varies more than 10% when using either of 



The species most sensi- 
24 Ms 



the rates from [68(, [42J, or |61|. 

tive to the changes in the 20 Ne + a ^ z ^Mg + 7 rates is 
24 Mg, whose yield changes by 156% when using the rate 
from [42| instead that from jgl], and by 84% when using 
the rate from [6l|. Several other species, like 27 Al and 
30 Si, experience changes in the range 70-90%. The yields 
obtained using the theoretical rate in (57j (not shown in 
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Table IXXIII) agree quite well with those obtained using 
the experimental rate in 61| . On the contrary, using the 
theoretical rates in [HI, obtained with the same code as 
in (57j but with different nuclear inputs, gives yields that 
differ from those belonging to the rates from [6l| by as 
much as 42%. 

The three most recent evaluations of the rate of the 
reaction 12 C + a 16 + 7 in JINA are from |6lj, [64|. 
and [HI . We have recomputed the supernova nucleosyn- 
thesis for these three prescriptions of the rate. Not a 
single species experiences a change of abundance larger 
than 10%. 

The reaction 58 Ni + a <=± 62 Zn + 7 is important for 
the alpha-rich freeze-out of NSE, which affects a large 
portion of SNIa ejecta. Hence, we have recomputed the 
nucleosynthesis with the three most recent rates given 
in JINA, from |42| and from (5?J using different nuclear 
inputs. Not a single species experiences a change of abun- 
dance larger than 10%. 



5. (a, n) reactions 

The different prescriptions for the reactions that have 
any \Di\ > 0.3 in Table IXl give rates that agree with 
each other within the factor often explored in this paper. 
These reactions are: 17 + a 20 Ne + n, 21 Ne + a «=i 
24 Mg + n, 25 Mg + a <=t 28 Si + n, 26 Mg + a <=t 29 Si + n, 
27 Al+a <=» 30 P+n, 29 Si+a «=i 32 S+n, 30 Si+a <=* 33 S+n, 
and 34 S + a T± 37 At + n. The rates are from a variety of 
sources, most of them from calculations with the NON- 
SMOKER code using different nuclear inputs, but there 
are as well rates based on experimental measurements by 

El El El. 



, p) reactions 



The different prescriptions for the rates of the reactions 



17 Ne + a ; 
26 Mg + p 



: 20 Na 
. . 27 Al ■+ 
Cr + a *± 



P- 



20 



Ne- 

30 c 



- Q 



23 Na + p, ^ 3 Na + a <=* 
a f= -Si + p, 34 S + a *± 37 C1 + p, 
50 Mn + p all agree within a factor of 



and 

ten for the temperature range of interest to us. These 
reaction rates come from several calculations with the 
NON-SMOKER code using different nuclear inputs, as 
well as experimental measurements in [6ll EH E3 ■ The 
different rates of the reactions 58 Ni + a <=^ 61 Cu + p and 



Ni- 
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Cu + p, both of importance for the alpha- 



rich freeze-out from NSE, agree quite well within the tem- 
perature range of interest. Thus, we do not continue with 
the analysis of these reaction rates. 

The three most recent evaluations of the rate of the 
reaction 24 Mg + a <=^ 27 A1 + p in the JINA database are 
from [H, EH E3 ■ Within the temperature range of inter- 
est, the rate from the last reference differs from the other 
rates by as much as three orders of magnitude. Hence, 
we have recomputed the nucleosynthesis for the three dif- 
ferent prescriptions of this rate. There is no significant 



difference in the chemical composition obtained with the 
rates of [42| and i.e. not a single species experi- 
ences a change of abundance larger than 10%, which is a 
consequence of the match between the rates from these 
two references. However, when comparing the yields ob- 
tained with the rate from [62j with those obtained with 
the rate from (60l | , many important variations in the yield 
of species with a significant abundance show up. 

The three most recent evaluations of the rate of the 
reaction 28 Si -I- a 31 P + p in JINA are from the same 
references as the reaction on 24 Mg. Although the discrep- 
ancies of the rates are not as large as for Mg, they reach 
one order of magnitude. Hence, we have recomputed as 
well the supernova nucleosynthesis for the three prescrip- 
tions of the rate of the reaction on 28 Si. However, not a 
single species experiences a change of abundance larger 
than 10%. 

The reaction 44 Ti + a 47 V + p is important for 
bridging the gap between QSE groups in silicon burn- 
ing. Hence, we have recomputed the nucleosynthesis with 
the three most recent rates given in JINA, from [42[ and 
from [57j using different nuclear inputs. Even though 
these rates differ up to near an order of magnitude in 
the range of temperatures of interest, not a single species 
experiences a change of abundance larger than 10%. 



7. Summary of the sensitivity of the yields to different rate 
prescriptions 



The most notable reactions in Table IXXIII are 30 Si + 
p i=> 31 p + 7j 20 Ne + a <=2 24 Mg + 7, and 24 Mg + a <=i 
27 A1 + p. Within the set of most influential reactions 
and most influenced species, there are two important de- 
tails. First, very few reactions appear whose parent nu- 
clei belong to the Fe-group or that are important for the 
bridging of the QSE groups in silicon burning. Second, 
in Table IXXIII there arc no product species belonging to 
the Fe-group whose yield depends significantly on the ex- 
plored reaction rates. This is most remarkable because 
the elements of the Fe-group constitute the main nucle- 
osynthetic products of SNIa. 



D. Sensitivity to different temperature ranges 

As a final step in our present study, we now perform an 
analysis of the temperature dependence of the sensitivity 
of the yields to the reaction rates. We analyze here the 
three most notable rates found in Section IIV C 71 To 
this end, we consider again a fixed enhancement factor 
of the rates, a factor of ten, but this time we limit it to 
a temperature window 10 9 K wide. We explore windows 
centered on temperatures 1.5, 2.5, 3.5, 4.5, 5.5, and 6.5 x 
10 9 K. We show the results in Figs. [T3] to [T5J where we 
plot for selected species the relative variation of their 
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yield, defined as 



(6) 



where m, is the mass of nucleus i ejected in our reference 
model (Table |n|, m\ is the yield of the same species when 
the rate of the reaction being analyzed is multiplied by a 
factor of ten independent of temperature, and my is the 
yield when the rate of the reaction is multiplied by the 
same factor only in the window j x 10 9 < T < (j + 1) x 
10 9 K. We selected the species to plot in the figures from 
among those with non-negligible abundances that present 
a large difference between rrt; and m^. We also required 
that the species covered a wide range of Z. Finally, we 
choose the same species to explore the sensitivities of all 
three reactions being considered: 20 Ne, 24 Mg, 26 Al, 30 Si, 
32 P, 35 S, 38 Ar, and 47 Ti. 
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FIG. 13. Sensitivity of the yield of selected species to varia- 
tions by a factor of xlO in the rate of the reactions 30 Si + p ±± 
P + 7 in different temperature ranges of size 10 9 K. The 
points give the difference between the yield of a species for 
an enhanced rate in a temperature window and its yield in 
our reference model, normalized by the difference between the 
yield for an enhanced rate at all temperatures and the yield 
of our reference model, see Eq. [6] The points are centered 
on each temperature window, and each symbol represents a 
product nucleus as follows: open triangles stand for 20 Ne, 
solid triangles for 24 Mg, crosses for 26 Al, open pentagons for 
30 Si, solid pentagons for 32 P, stars for 35 S, asterisks with seven 
vertices for 38 Ar, and open circles for 47 Ti. The dashed lines 
(scaled according to the right axis) give the logarithm of the 
ratio of the different rates of the reaction 30 Si + p — > 31 P + 7 in 
JINA. The short-dashed line belongs to the ratio of the rate 
from [42l ] to that from [fjtj , while the long-dashed line belongs 
to the ratio of the rate from to that from [gj}. The hor- 
izontal solid line marks the zero of axes, i.e. no variation of 
the yield and rate ratio equal to one. 



We show in Fig. [13] the results for the reaction 
30 Si + p <=* 31 P + 7. We note that all the nuclei follow 
the same behavior with respect to the temperature win- 
dow in which the reaction rate is modified. There is a 
modest variation of order 20% of the yields for the win- 
dow centered on 1.5 x 10 9 K, and a large increase for the 
next window, centered on 2.5 x 10 9 K. In the window 
centered on 3.5 x 10 9 K there are some species which ex- 
perience large variations of their yields while others are 
scarcely affected at all. The yield of 20 Ne (open trian- 
gles) shows a peculiar behavior, with nj < in the first 
thermal window, meaning that increasing the rate of the 
reaction only at low temperatures (T < 2 x 10 9 K) re- 
sults in a variation of the yield of 20 Ne of opposite sign 
as that obtained if the reaction rate is increased for any 
temperature. Note that the sign of of 20 Ne in the next 
window is positive, and it has the largest r»j among the 
species shown in the figure: increasing the reaction rate 
only in the interval 2 x 10 9 < T < 3 x 10 9 K produces 
a change of the yield of this species that is as much as 
that obtained by increasing the reaction rate for all tem- 
peratures. Modifying the rate on thermal windows above 
4 x 10 9 K has no effect on any of the final abundances of 
the species. 

We show as well in Fig. Q~3] the ratio of the rates be- 
longing to the three prescriptions adopted for the rate of 
the reaction 30 Si + p ^ 31 P + 7, which were discussed 
in Section IIV C 21 and in Table IXXIII The uncertainty 
in the rates derived from these different prescriptions is 
more or less uniform for temperatures above ~ 2x 10 9 K. 
The rate from [57J differs most from that based on [60] at 
temperatures where the yields are most sensitive to this 
reaction rate. However, the discrepancy between these 
rates is much lower than the factor of ten used in our 
simulations, thus we believe that the supernova yields 
should not be affected by any reasonable future change 
of this reaction rate. 

Figure [14] summarizes the results for the reaction 
20 Ne + a ^ 24 Mg + 7. The most noticeable difference 
with respect to Fig. [T3] is the behavior and range of the 
variations of the yield of 38 Ar (asterisks). The maximum 
sensitivity of this species occurs in the temperature win- 
dow 4 x 10 9 < T < 5 x 10 9 K, where the change of its yield 
reaches a value seven times larger than the change with 
a rate modified at all temperatures. This is compensated 
by the fact that modifying the rate at temperatures in 
the interval 2 x 10 9 < T < 4 x 10 9 K produces a change 
of the yield of 38 Ar of opposite sign. The rest of nuclei 
plotted show a behavior similar to the one in Fig. [T3l 
with maximum |ry| ~ 1.5 ( 32 P, solid pentagons). 

As revealed by Fig. Q3J the different prescriptions for 
the rate of the reaction 20 Ne + a — > 
maximum discrepancy by a factor of * 
ature range 10 9 - 10 10 K. However, both the rate from 
[42| and that from [6l| differ from the rate given by [68| 
by a similar factor in the interval 2xl0 9 <T<4xl0 9 K. 

Finally, we show in Fig. [15] the results for the reac- 
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FIG. 14. Sensitivity of the yield of selected species to varia- 
tions by a factor of x 10 in the rate of the reactions 20 Ne+a «=* 
24 Mg + 7 in different temperature windows of size 10 9 K. 
The meaning of the points is the same as in Fig. 1131 The 
dashed lines (scaled according to the right axis) give the 
logarithm of the ratio of the different rates of the reaction 
20 Ne + a -s> 24 Mg + 7 in JINA. The short-dashed line belongs 
to the ratio of the rate from (4^] to that from [6?§], while the 
long-dashed line belongs to the ratio of the rate from [6lJ to 
that from [68j]. The two horizontal lines mark the zero of the 
left axis, i.e. no variation of the yield (solid line), and the zero 
of the right axis, i.e. rate ratio equal to one (dotted line). 



FIG. 15. Sensitivity of the yield of selected species to 
variations by a factor of xlO in the rate of the reactions 
24 Mg + a <=i 27 Al + p in different temperature windows of size 
10 9 K. The meaning of the points is the same as in Fig. 1131 
The dashed lines (scaled according to the right axis) give the 
logarithm of the ratio of the different rates of the reaction 
24 Mg + a -> 27 Al + p in JINA. The long-dashed line belongs 
to the ratio of the rate from [Z^] to that from [6(| , while the 
short-dashed line belongs to the ratio of the rate from [62] to 
that from The two horizontal lines mark the zero of the 
left axis, i.e. no variation of the yield (solid line), and the zero 
of the right axis, i.e. rate ratio equal to one (dotted line). 



of 30 Si (open pentagons) whose yield experiences varia- 
tions up to seven times larger in the temperature win- 
dow 3 x 10 9 < T < 4 x 10 9 K than when the rate is 
modified for all temperatures. Furthermore, the vari- 
ation of its yield changes sign if the thermal window 
is 2 x 10 9 < T < 3 x 10 9 K, still reaching \nj\ ~ 7. 
The rest of the nuclei show a behavior similar to that in 
Fig-HSl although their maximum |r»j | is now a bit larger, 
max (\nj |) ~ 2. 

The ratio of the rates of the reaction 
24 Mg + a <=t 27 A1 + p from [42] and agree quite 

well for T > 2 x 10 9 K (see Fig. Q~5])'. On the other 
hand, the rate from [(52[ differs from the other two by 
more than two orders of magnitude at high tempera- 
tures. However, in the range of temperatures where the 
abundances plotted are most sensitive to this reaction, 
2xl0 9 <T<4xl0 9 K, their discrepancy is less than 
a factor of ~ 30. It is interesting to note that, using 
these rates, the change of the sensitivity of the yield 
of 30 Si in the two temperature windows with |ry| ~ 7 
almost compensates each other, with the result that the 
final yield of this species is negligibly affected by using 
the rate from [62| instead of that from [60| (thus, it 
does not appear in the row reserved for this reaction in 



Table IXXHi- 



V. CONCLUSIONS 

We have computed the chemical composition of a refer- 
ence SNIa model with a nucleosynthetic post-processing 
code that takes as inputs the nuclear data and the ther- 
modynamic trajectories of each mass shell. Our reference 
SNIa model is the one-dimensional delayed-detonation 
model DDTc in |48[ . Our nucleosynthetic calculations in- 
clude 3138 reactions during the integration of the nuclear 
evolutionary equations but only 1096 of them can con- 
tribute significantly to the nucleosynthesis of the super- 
nova model. In this paper, we have explored the sensitiv- 
ity of the SNIa explosive nucleosynthesis to simple varia- 
tions on nuclear reaction rates (either a fixed enhance- 
ment factor or one that decreases monotonously with 
temperature) and comparisons between different theoret- 
ical and experimental prescriptions of the rates. 

The nucleosynthesis resulting from our Type la super- 
nova model is quite robust with respect to variations of 
nuclear reaction rates, with the exception of the fusion of 
two 12 C nuclei. The energy of the explosion changes by 
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less than ~ 4% when the rates of the reactions 12 C + 12 C 
or 16 + 16 are multiplied by a factor of 10 or 0.1. The 
changes in the nucleosynthesis due to the modification 
of the rates of these fusion reactions are also quite mod- 
est, for instance no species with a mass fraction larger 
than 0.02 experiences a variation of its yield larger than 
a factor of two. The robustness of the production of 56 Ni 
and many other Fe-group isotopes stands out. If the en- 
hancement factor of the rates is a decreasing function of 
temperature, the effect on the yields of all the species 
is even less than with a fixed enhancement factor. For 
instance, no species experiences a variation of its yield 
larger than 30% with respect to the rate of 16 + 16 
or with respect to the 3a rate. We have checked that 
the modifications in the nucleosynthesis produced when 
the fusion rates are modified only in the nucleosynthetic 
code are quite similar to those obtained when the rates 
are modified in the full supernova simulation. 

We provide the sensitivities of the yield of each rele- 
vant nuclear species ejected in the supernova with respect 
to those nuclear reactions that affect it. In general, the 
yields of Fe-group nuclei are less sensitive than the yields 
of intermediate-mass elements. However, the yields of 
28 Si, and 32 S, as well as 54 Fe, 56 Ni, and 58 Ni do not 
change appreciably within the range of enhancement fac- 
tors of the nuclear reaction rates explored here. The only 
reactions for which the relative change of the abundance 
of any species is larger (in absolute value) than the rela- 
tive change in the rate (|A| > 1) is 30 Si + p ^ 31 P + 7. 
In general, radiative captures of protons are the group 
of reactions with the largest influence on the supernova 
yields. Other important groups of reactions are radia- 
tive captures of a particles, most notably the reactions 
20 Ne + a 24 Mg + 7, for which there are 33 species 
whose yields change by more than 12%, and the (a, p) 
reactions 13 N + a <=* 16 + p, 20 Ne + a ^ 23 Na + p, 
23 Na + a <=* 26 Mg + p, 24 Mg + a z± 27 A1 + p, and 
27 A 1 + a «=♦ 30 S i + p. 

We have discussed the sensitivity of the nucleosynthe- 
sis to the rates of reactions that take part in the most 
relevant nucleosynthetic processes in SNIa. Modifying 
the rates of the reactions that bridge the gap between 
QSE groups in explosive silicon burning has a very lim- 
ited effect on the yields. They affect mainly nuclei that 
belong to this gap. Changing the rate of reactions rele- 
vant for the alpha-rich freeze-out does not produce im- 
portant changes on the abundances of nuclear species ei- 
ther. This can be explained both by the small amount 
of matter that goes through alpha-rich freeze-out from 
NSE and by the small excess of alpha particles at freeze- 
out (see Fig. [T] where the maximum value of Y a does 
not attain 0.01 mol g _1 ). It is as well remarkable that 
the reactions involving nuclei with Z > 22 have a tiny 
influence on the supernova nucleosynthesis. 

We have relied on the JINA REACLIB Database to 
estimate realistic uncertainties of the most relevant re- 
action rates, by comparing the most recent prescrip- 
tions for the rates of these reactions. We have paid 



special attention to the reactions 30 Si + p <^ 31 P + 7, 
20 Ne + a <=± 24 Mg + 7, and 24 Mg + a £± 27 Al + p, espe- 
cially the last one for which there is a discrepancy of up 
to three orders of magnitude between the rates due to 
[60| and (62J. In spite of this large difference of rates, the 
maximum change in the yields is only 52%, belonging to 
35 S. 

There are two main reasons for the small relative im- 
pact of the uncertainties of individual nuclear reaction 
rates on the supernova yields. First, the nuclear flows 
that determine the final abundances during the super- 
nova explosion are driven collectively by many reactions, 
which are much faster than the hydrodynamic explosion 
timescale because of the high temperatures involved. The 
relevance of any individual rate is much diluted within 
this large pool of reactions. A similar conclusion was 
reached by [H| in the context of Type II supernovae. 
They cite three major causes, which we can adapt to nu- 
cleosynthesis in SNIa: 1) the dominant nuclear flows are 
governed by the fusion reactions of the fuel, carbon and 
oxygen, while the rest are only perturbations on the main 
stream, 2) the nuclear flow follows the path of least re- 
sistance, i.e. if one reaction rate drops by a large factor 
there is always another reaction capable of playing its 
role, and 3) if the freeze-out from high temperatures is 
fast enough, the rates of individual reactions are much 
less important than the properties of nuclei (binding en- 
ergy, partition function). 

Second, there are narrow temperature ranges where 
the yields are more sensitive to the rates. For instance, 
the temperatures at which a modification of the rate of 
the above-mentioned three reactions has a larger impact 
are in the range 2 x 10 9 < T < 4 x 10 9 K (see Figs. [TBI to 
fTS")) . One kind of rate uncertainty we have not explored is 
that due to the erroneous location of a resonance. Such 
a kind of error might originate an increase of the rate 
(with respect to the presently recommended one) in a 
temperature range and a decrease in a contiguous one. 
If this were the case, the changes of the yields of some 
species might be exacerbated. Thus, this kind of error 
in the nuclear reaction rates might be the most relevant 
with respect to the supernova yields. 

We conclude that the explosion model chiefly deter- 
mines the element production of Type la supernovae, 
and derived quantities like their luminosity, while the in- 
dividual nuclear reaction rates used in the simulations 
have a small influence on the kinetic energy and final 
chemical composition of the ejecta. Often, it is argued 
that discrepancies of up to a factor of two between iso- 
topic ratios in SNIa ejecta and those in the solar sys- 
tem, especially within the Fe-group, can be attributed 
to uncertainties in nuclear reaction rates. Our results 
show that the uncertainty in individual thermonuclear 
reaction rates cannot account for this factor. It remains 
to be seen if the yields are more sensitive to uncertain- 
ties in nuclear masses, weak interaction rates, or to the 
simultaneous modification of the bulk of thermonuclear 
reaction rates. The sensitivity of the supernova nuclc- 
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osynthesis to simultaneous random modifications in the 
bulk of thermonuclear reaction rates will be the subject 
of future work. In this respect, it is interesting our find- 
ing that the most influential reactions depict a clear path 
in a plot Z vs. A (Fig. IT21) . going from 12 C up to 37 Ar 
through many branches involving mainly reactions with 
a particles plus the fusion reaction 12 C + 12 C. Modi- 
fications of these rates 'in phase' may have interesting 
consequences for the chemical composition of supernova 
ejecta. 

Finally, it is worth noting that reaction rate variations 
may also have an impact on the hydrostatic evolution of 
the progenitor of the exploding white dwarf. Given the 
robustness of the explosive yields, it may well be that 
changes in progenitor evolution are the largest source of 



reaction rate sensitivity in thermonuclear supernovae. 
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